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PREFACE 

One  of  the  most  pressing  problems  of  Air  Force  Installation  planners  in  the 
1970's  is  minimizing  environmental  pollution.  Air  Weather  Service  personnel 
can  give  valuable  assistance  to  these  planners  by  providing  realistic  estimates 
of  atmospheric  diffusion.  This  Technical  Report  has  been  prepared  as  an  aid  to 
making  diffusion  estimates.  The  scope  of  this  report  is  limited  to  calculating 
pollution  concentrations  from  on-base  sources. 
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Chief  Scientist 

USAF  Environmental  Technical  Applications  Center 
Washington,  D.  C. 
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Chapter  1 
INTRODUCTION 


1.  General. 

Air  pollution  is  not  a  new  phenomenon  and  has  been  the  concern  of  many  nota¬ 
ble  Individuals  and  scientists  for  years.  Nevertheless,  air  pollution,  along 
with  water  pollution,  has  suddenly  become  a  subject  of  national  importance.  In 
the  last  few  years  several  Presidential  Executive  Orders  have  been  issued  with 
the  aim  of  reducing  sources  of  air  pollution  at  Federal  installations.  Thus, 
the  Air  Force  is  deeply  concerned  with  reducing  air  pollution  around  their  in¬ 
stallations  and  has  embodied  this  concern  in  AFR  161-22.  Provisions  of  this 
regulation  task  Air  Weather  Service  with  providing  meteorological  support  to 
air-pollution  oontrol  and  abatement  personnel.  As  a  consequence,  detachment 
commanders  and  staff  weather  officers  can  expect  an  increasing  number  of  in¬ 
quiries  from  civil  engineering  and  environmental  health  personnel  about  the  me¬ 
teorological.  aspects  of  air  pollution.  Intense  Interest  in  air  pollution,  how¬ 
ever,  is  not  limited  to  these  two  activities  only.  All  aspects  of  operations 
and  training  within  the  Air  Force  must  operate  with  an  awareness  of  the  pollu¬ 
tion  problem.  For  example.  Headquarters  USAF  has  recently  instructed  the  major 
air  commands  to  conduct  fire-fighting  training  during  periods  of  good  atmos¬ 
pheric  dispersion. 

Historically,  military  interests  in  atmospheric  dispersion  dates  back  to  the 
use  of  toxic  gases  in  World  War  I.  Much  of  the  early  work  on  atmospheric  dis¬ 
persion  was  conducted  by  British  scientists  in  connection  with  chemical  warfare 
research.  Later  American  researchers  also  recognized  the  importance  of  such 
investigations.  The  emphasis  placed  by  Air  Weather  Service  on  atmospheric  dis¬ 
persion  began  in  earnest  with  the  recognition  that  liquid  fuels  for  Titan  mis¬ 
siles  were  toxic  and  that  a  potential  safety  hazard  existed  with  respect  to 
storage  and  transfer  of  these  potentially  dangerous  fuels.  The  assistance  re¬ 
quired  of  Air  Weather  Service  personnel  by  the  operational  aspects  of  the  Titan 
missile  complexes  intensified  the  research  for  knowledge  of  dispersion  within 
the  free  atmosphere  next  to  the  earth's  surface. 

In  order  to  ably  fulfill  their  new  responsibilities  as  advisors  to  the  anti- 
pollution  effort,  AWS  forecasters,  and  especially  DETCOs  and  STAFFMETS,  should 
become  familiar  with  the  natural  processes  that  control  the  dispersion  of  con¬ 
taminants.  This  technical  report  is  intended  to  give  the  field  forecaster  an 
appreciation  of  atmospheric  dispersion  processes  and  practical  methods  for 
assessing  the  ability  of  the  air  to  dilute  and  carry  away  contaminants.  A  com¬ 
prehension  of  such  processes  entails  a  knowledge  of  the  very  lower  layers  of 
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the  atmosphere. 

2 »  Various  Scales  of  Air  Pollution. 

There  are  three  broad  scales  of  air  pollution  which  can  be  defined  on  the 
basis  of  geographical  extent  and  degree  of  dependence  on  local  meteorological 
conditions.  In  order  of  decreasing  geographical  scale,  they  are  global  pollu¬ 
tion,  area  pollution,  and  local  pollution. 

a.  Qlobal  Air  Pollution.  On  this  scale,  while  local  effects  are  of  little 
importance,  the  general  circulation  imposes  some  restriction  on  the  ultimate 
spread  of  pollutants.  The  principal  factors  controlling  the  concentration  of 
pollutants  on  this  scale  are  the  rate  at  which  the  pollutants  are  being  intro¬ 
duced  into  the  atmosphere  worldwide  and  the  "residence  time,"  i.e.,  the  length 
of  time  it  takes  for  a  pollutant  to  settle  out  as  particles,  be  washed  out  by 
precipitation,  to  diffuse  upward  and  eventually  into  space,  or  to  be  converted 
through  some  chemical  process  into  an  invisible,  harmless  gas. 

b.  Area  Pollution.  A  second  scale  of  air  pollution  is  typified  by  the  Los 
Angeles  "smog"  problem.  Here,  existing  meteorological  conditions  combined  with 
unique  terrain  features  act  to  confine  the  pollutants  in  a  localized  region. 

The  meteorological  parameters  of  importance  on  this  scale  are  dealt  with  rou¬ 
tinely  in  forecasting,  viz.,  wind  speed  and  temperature  inversions.  Again,  as 
on  the  global  scale,  the  details  of  diffusion  from  individual  sources  are  not 
too  important.  This  is  the  type  and  scale  of  air  pollution  considered  in  the 
air  pollution  potential  forecasts  that  are  relayed  to  AWS  detachments  in  the 
continental  United  States  on  the  communication  circuits. 

c.  Local  Pollution.  On  the  smallest  scale  the  concentration  of  pollutants 
mostly  depends  on  the  details  of  the  processes  of  atmospheric  diffusion.  An 
example  of  this  type  of  pollution  is  the  zone  of  acrid  fumes  that  often  exists 
downwind  from  pulp  mills.  This  type  of  local  phenomenon  could  also  be  contrib¬ 
uting  to  both  of  the  other  scales  of  air  pollution  but  at  levels  that  would  not 
be  readily  noticed.  An  air  pollution  source  on  the  local  scale  is  usually  con¬ 
sidered  serious  only  if  it  produces  concentrations  of  contaminants  that  exceed 
certain  limits  within  populated  areas.  However,  as  will  be  explained  in  the 
next  section,  the  limits  on  concentration  apply  to  the  total  concentration, 
i.e.,  the  sum  of  all  local  and  area  sources. 

3.  Air  Pollution  Limits  and  the  Concept  of  Air  Quality. 

Most  of  the  state  and  local  air  pollution  laws  express  their  limits  in  terms 
of  parts  of  contaminant  per  million  parts  of  contaminated  air  at  1013-mb  pres¬ 
sure  and  25° C.  The  general  concept  of  enforcement  is  to  evaluate  the  concentra¬ 
tion  at  the  downwind  property  line.  In  general,  this  concentration  will  be  the 
sum  of  the  local,  pollution  introduced  at  the  installation  (diluted  in  the  dis¬ 
tance  from  release  point  to  the  property  line),  plus  whatever  pollution  is  pres¬ 
ent  in  the  area  due  to  other  sources.  This  residual  pollution  from  other 
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eouroee  ia  known  as  the  background  level. 

It  should  be  apparent  now  that  the  background  pollution  level  ia  aa  impor¬ 
tant  as  the  local  diffusion  characteristics  in  determining  and  controlling  the 
amount  of  pollutants  that  can  safely  be  emitted.  When  the  background  is  high, 
only  a  small  amount  of  pollutant  can  be  released  before  a  critical  condition  is  . 
reached,  even  if  diffusion  conditions  are  locally  good. 

Air-quality  standards  fall  into  two  general  types  -  long  term  and  short  term. 
The  long-term  limits  apply  to  average  concentrations  measured  over  periods  of 
several  hours  to  a  month.  Short-term  limits,  which  apply  to  situations  lasting 
an  hour  or  less,  have  not  yet  been  adopted.  However,  the  Environmental  Protec¬ 
tion  Agency  will  begin  issuing  short-term  guidelines  in  1971.  State  and  local 
agencies  can  be  expected  to  adqpt  short-term  limits  very  soon  thereafter. 

Air  Weather  Service  personnel  normally  would  not  become  involved  in  setting 
air-quality  standards.  However,  it  is  within  the  framework  of  these  standards 
or  limits  that  the  staff  meteorologist  must  make  his  pollution  forecasts.  The 
situation  is  very  similar  to  setting  critical  crosswind  components  for  aircraft 
landings.  Two  similarities  are  apparent.  First,  the  limits  are  different  for 
each  pollutant  Just  as  different  crosswind  components  apply  to  each  aircraft 
type.  Second,  the  setting  of  the  pollution  limits  is  a  subjective  process  in 
which  medical  or  industrial  hygiene  personnel  evaluate  past  data  on  exposures  to 
pollutants  Just  as  flying  safety  personnel  use  their  Judgment  and  the  accumu¬ 
lated  data  on  accidents  to  set  critical  crosswinds.  In  both  cases,  the  opera¬ 
tional  impact  may  be  considered  in  balancing  the  requirements  of  health  and 
safety  against  the  requirements  of  mission  accomplishment.  However,  in  the  case 
of  air  pollution,  the  balance  must  be  made  at  a  higher  level  of  government  since 
the  health  and  safety  of  the  general  public,  not  Just  military  personnel  and 
equipment,  are  involved. 

4.  The  Relationship  Between  Meteorology  and  the  Concentration  of  Air  Pollutants. 

Meteorology  and  air  pollution  have  a  direct  relationship.  The  occurrence  of 
air  pollution  in  any  degree  depends  on  the  presence  or  absence  of  dilution  and 
diffusion  processes.  In  the  free  atmosphere,  both  of  these  processes  depend  al¬ 
most  entirely  on  familiar  meteorological  parameters,  viz.,  turbulence,  wind, 
temperature,  precipitation,  atmospheric  stability,  inversions,  etc.  This  rela¬ 
tionship  can  best  be  illustrated  by  two  examples.  For  instance,  consider  the 
warning  found  on  the  containers  of  many  kinds  of  cleaning  fluids  or  solvents: 
"Caution  —  Avoid  excessive  breathing  of  fumes.  Use  only  in  a  well -ventilated 
place."  Why  is  the  ventilated  place  less  dangerous  than  a  closed  room?  The 
reason  is  the  presence  of  the  dilution  process,  yhe  ventilated  area  allows 
fresh  air  to  mix  with  and  dilute  the  toxic  fumes  of  the  cleaning  fluid  so  that 
the  remaining  concentration  is  not  harmful.  The  property  of  dilution  is  eco¬ 
nomically  one  of  the  most  useful  properties  of  the  atmosphere.  In  our  daily 
living  countless  tons  of  waste  products  annually  are  disposed  of  simply  by 
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discharging  the  material  into  the  atmosphere  and  letting  the  winds  carry  it 
away.  As  the  waste  material  is  transported,  it  la  also  mixed  through  larger  and 
larger  volumes  of  air  by  the  action  of  turbulent  eddies  with  the  concentration 
beooming  weaker  and  weaker.  This  is  an  example  of  the  way  in  which  the  natural 
environment  makes  the  handling  of  toxic  chemicals  safer. 

A  second  example,  however,  shows  how  these  same  atmospheric  processes,  pre¬ 
viously  shown  to  be  so  advantageous,  can  make  other  situations  more  dangerous. 
Suppose  a  railroad  tank-car  full  of  chlorine  is  derailed  and  ruptures.  Police 

and  Civil  Defense  officials  rush  to  evacuate  people  from  the  area  immediately 

downwind  from  the  wreck.  Why?  Because  the  wind  will  carry  the  deadly  chlorine 
fumes  from  the  spot  of  the  accident  to  other  places  as  much  as  several  miles 
away.  With  insufficient  diffusion  this  brings  danger  where  none  previously  ex¬ 
isted.  Suppose  the  air  were  absolutely  still.  With  no  wind  there  would  be  no 
spread  of  dangerous  fumes  except  by  the  comparatively  slow  process  of  molecular 

diffusion  and  thus  no  other  area  would  be  dangerously  contaminated  by  the  chlor¬ 

ine. 

In  every  real  situation,  the  good  and  the  bad  are  both  represented.  For,  as 

in  the  example  of  the  last  paragraph,  the  same  moving  atmosphere  that  created  a 

danger  by  transporting  the  chlorine  also  helps  reduce  the  danger  by  diluting 
the  concentration  of  the  gas  as  it  moves.  The  essence  of  the  relationship  be¬ 
tween  the  atmosphere  and  the  release  of  toxic  substances  can  be  found  in  the 
following  question.  Will  the  mixing  action  of  the  atmosphere  dilute  the  mate¬ 
rial  to  a  tolerable  level  fast  enough  to  keep  the  spread  from  being  a  danger? 

This  question  is  indeed  a  challenging  one  because  the  cleansing  efficiency  of 

the  atmosphere  varies  from  place-to-place,  day-to-dsy,  or  even  hour-to-hour. 

The  same  principles  apply  to  nontoxic  air  pollution.  That  is,  the  winds 
carry  pollution  from  the  source  to  other  places  while  diluting  it  along  the  way. 
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Chapter  2 

FUNDAMENTALS  OF  MCQKOMETEOROLOGY 

I  1  !(  . 


1.  Qenaral. 

The  term  micrometeorology  Is  used  In  this  report  in  Its  broadest  sense]  that 
is,  as  the  meteorology  of  the  layer  from  the  earth's  surface  to  an  altitude 
where  the  effects  of  the  immediate  underlying  surface  on  the  temperature  and 
motion  of  the  atmosphere  become  negligible,.  This  definition  of  micrometeorology 
includes  a  considerable  portion  of  what  is  often  considered  as  mesometeorology. 

Most  Air  Weather  Service  meteorologist o  have  received  the  bulk  of  their  me¬ 
teorology  experience  as  synoptic  meteorologists.  As  such,  they  were  trained  to 
ignore  or  "filter  out"  of  their  analyses  any  local  irregularities  that  existed 
on  a  scale  of  less  than  50  miles  or  so,  except  for  incorporating  local  effects 
in  terminal  forecasts.  This  "filtering"  is  quite  appropriate  for  depicting  the 
general  state  of  the  weather  and  for  predicting  those  aspects  of  weather  that 
change  according  to  the  movement  of  major  pressure  systems.  On  the  other  hand, 
the  air  motions  that  affect  the  transport  and  diffusion  of  air  pollutants  are 
much  more  local  and  depend  on  atmospheric  conditions  a  few  hundred  feet  in  ver¬ 
tical  extent  and  hundreds  of  yards  to  a  few  miles  in  the  horizontal.  When  one 
is  considering  the  transport  and  dilution  of  pollutants,  the  most  important 
variables  are  wind  direction  and  temperature  lapse  rate.  In  that  near-ground 
layer  both  the  temperature  lapse  rate  and  the  air  movements  exhibit  character¬ 
istic  variations  from  day  to  night  and  from  season  to  season.  Furthermore, 
small-scale  differences  in  the  temperature  and  wind  fields  can  be  attributed  to 
local  irregularities  in  the  underlying  terrain.  This  chapter  will  review  the 
physical  basis  for  the  normal  variation  of  wind  and  temperature  near  the  ground. 

2.  The  Physical  Basis  of  Meteorological  Variation  Near  the  Ground. 

Over  level  terrain,  both  the  normal  diurnal  variation  and  the  local  irregu¬ 
larities  of  the  temperature  and  wind  fields  can  be  explained  primarily  in  terms 
of  the  response  of  the  underlying  surface  to  incoming  solar  radiation  by  day  and 
outgoing  Infrared  radiation  at  night. 

a.  The  Balance  of  Energy  at  the  Surface.  The  balance  of  energy  at  the 
earth's  surface  in  the  daytime  is  roughly  shown  in  Figure  la.  The  wide  arrow 
labeled  "net  Incoming  solar  radiation"  is  the  resultant  of  all  radiation  proc¬ 
esses  (absorption,  reflection,  reradiatlon)  acting  at  the  ground  surface.  About 
half  of  the  net  Incoming  radiant  energy  is  conducted  downward  into  the  ground 
during  the  morning  as  the  ground  warms.  This  point  will  be  elaborated  later. 

The  remaining  half  of  the  energy  is  split  between  evaporation  and  upward  heat 
transport.  The  relative  magnitudes  of  these  two  processes  depends  mostly  on  the 


Technical  Report  214  April  1971 

moisture  available  for  evaporation. 

The  energy  flow  at  the  surface  during  the  nighttime  is  depicted  in  Figure  lb. 
At  night,  heat  flows  from  the  ground  to  supply  much  of  the  energy  that  la  lost 
from  the  ground  surface  as  effective  outgoing  radiation.  The  relative  magnitudes 
of  the  heat  flows  are  roughly  indicated  by  the  width  of  the  arrows,  which  are 
also  scaled  to  show  magnitudes  as  compared  with  Figure  la.  The  main  difference 
between  day  and  night  energy  transfer  is  that  heat  flows  into  the  ground  and  in¬ 
to  the  lowest  layers  of  the  air  by  day  and  the  flow  is  reversed  by  night.  These 
flow  directions  result  in  substantial  changes  from  day  to  night  in  the  tempera¬ 
ture  of  the  ground  surface  and  the  air  layer  nearest  it.  The  amount  of  temper¬ 
ature  change  depends  on  the  physical  characteristics  of  the  underlying  surface 
and  on  the  operation  of  all  three  modes  of  heat  transfer  within  the  atmosphere. 


NIT  IN60NIMI  SOLAN 
NADIATION 


Figure  1,  Energy  Balance  at  the  Ground 
(a)  Middle  of  the  Day,  (b)  At  Night. 


b.  Physical  Characteristics  Affecting  the  Response  of  the  Underlying 
Surface .  One  factor  affecting  the  energy  balance  at  the  ground  surface  is  the 
albedo,  or  reflectivity,  of  the  underlying  surface.  Those  surfaces  with  low  al¬ 
bedo  will  absorb  more  of  the  Incoming  solar  radiation.  For  example,  the  high 
midday  temperature  of  a  black-top  road  (low  albedo)  compared  to  a  nearby  field 
(higher  albedo)  is  primarily  due  to  this  albedo  difference.  However,  the  al¬ 
bedo  is  not  too  often  the  controlling  parameter.  Concrete  pavement  even  with 
its  higher  albedo  will  often  become  hotter  than  black-top.  In  this  case  the 
difference  is  found  primarily  in  the  heat  conductivity  and  the  heat  capacity 
(the  amount  of  heat  required  to  raise  the  temperature  of  a  substance)  per  unit 
volume  of  the  particular  surface.  The  heat  capacity  per  unit  volume  i3  the 
product  of  specific  heat  and  density.  If  the  density  and  specific  heat  of  the 
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underlying  surface  are  large,  a  great  amount  of  energy  must  be  absorbed  or  given 
off  at  the  surface  in  order  to  change  the  temperature  appreciably,  if  the  heat 
conductivity  of  the  surface  is  high,  the  heat  energy  will  flow  rapidly  into  or 
out  of  the  ground.  This  rapid  heat  flow  will  result  in  a  small  temperature 
change  occurring  over  a  large  depth,  while  a  low  heat  conductivity  would  result 
in  large  temperature  changes  through  a  shallow  depth. 

c.  Physical  Processes  Affecting  the  Response  of  the  Air  Layer  Near  the 
Ground.  We  learn  in  elementary  physics  courses  that  there  are  three  r  rocesses 
that  transfer  heat:  radiation,  conduction,  and  convection.  The  first  two  have 
already  been  mentioned  in  connection  with  the  variation  of  the  ground  tempera¬ 
ture.  Both  radiation  and  conduction  contribute  to  heat  transfer  in  the  atmos¬ 
phere,  but,  except  for  a  very  thin  layer  next  to  the  ground,  both  of  these  proc¬ 
esses  are  less  effective  for  heat  interchange  than  convection  in  its  broadest 
sense.  In  this  instance,  convection  means  transport  of  heat  due  to  any  motion 
of  fluid  particles  other  than  molecular.  Thus,  convection  includes  stirring  and 
mixing  by  all  causes  as  well  as  thermally-forced  convection,  which  we  associate 
with  dust  devils  or  cumulus  clouds.  The  total  heat  transferred  to  or  from  the 
ground  surface  by  all  of  these  processes  acting  in  the  air  is  less  than  the  heat 
transferred  by  any  one  of  the  various  processes  affecting  the  temperature  of  the 
ground.  Therefore,  heat  transfer  within  the  atmosphere  has  only  a  minor  influ¬ 
ence  on  the  range  of  temperature  at  the  ground  or  the  temperature  of  the  air 
immediately  next  to  the  ground.  Instead,  heat  transfer  within  the  atmosphere 
determines  how  deep  a  layer  will  respond  to  changes  of  ground  temperature  and 
the  speed  of  that  response.  As  an  extreme  case,  if  there  were  no  heat  transfer 
of  any  kind  in  the  air,  the  air  molecules  right  next  to  the  ground  would  acquire 
the  temperature  of  the  ground  and  the  temperature  of  the  rest  of  the  atmosphere 
would  riot  change.  In  reality,  there  is  considerable  heat  transport,  enough  that 
at  least  several  hundred  feet  of  the  atmosphere  experiences  a  very  noticeable 
change  of  temperature  from  day  to  night.  In  the  middle  of  the  day  and  late  at 
night,  when  the  air  and  ground  temperatures  are  changing  only  slowly,  the  tem¬ 
perature  gradients  will  be  inversely  proportional  to  the  effectiveness  of  the 
heat,  transfer. 

As  noted,  the  principal  mode  of  heat  transfer  in  the  atmospheric  surface 
layer  is  turbulent  mixing.  The  term  "turbulent,"  as  used  here,  refers  only  to 
the  randomness  of  the  motions  and  does  not  imply  the  violent  motion  that  a  fore¬ 
caster  normally  visualizes  when  he  considers  aircraft-type  turbulence.  Random 
motions  of  the  air  are  present  nearly  everywhere.  The  curling,  irregular  wisp 
of  smoke  from  a  cigarette  testifies  to  the  presence  of  some  turbulence  even  in¬ 
doors.  As  gentle  as  the  random  mixing  nay  be,  however,  it  is  still  10  to  100 
times  more  efficient  in  transferring  heat  than  is  molecular  conduction.  Out¬ 
doors,  where  the  mixing  is  more  vigorous,  turbulent  transfer  is  100  to  100,000 
times  more  efficient.  The  intensity  of  the  mixing  depends  on  the  speed  of  the 
wind,  the  roughness  of  the  terrain,  and  the  vertical  gradient  of  temperature. 
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The  effect  of  the  temperature  lapse  rate  depends  on  whether  It  Is  greater 
or  less  than  the  adiabatic  lapse  rate.  If  the  lapse  rate  Is  stable*  energy  will 
be  consumed  in  turbulent  mixing,  while  an  unstable  atmosphere  will  feed  energy 
Into  the  turbulent  motions.  This  characteristic  gives  the  atmosphere  a  pre¬ 
ferred  direction  of  heat  flow  -  upward.  If  the  atmosphere  io  extremely  un¬ 
stable,  the  intensity  of  turbulence  Is  great  and  a  large  amount  of  heat  Is  trans 
ferred,  which,  in  turn,  tends  to  reduce  the  lapse  rate.  On  the  other  hand,  the 
sharper  the  nocturnal  inversion  the  weaker  is  the  mixing  motion,  tending  to  in¬ 
tensify  the  inversion  when  radiation  processes  more  than  offset  the  weak  mixing. 

3*  Looal  Temperature  Structure. 

The  preceding  discussion  provides  a  basis  for  understanding  both  the  normal 
variations  of  temperature  near  the  earth's  surfaoe  and  the  differences  that 
arise  because  of  the  inhomogenieties  of  that  surface. 

a.  Diurnal  and  Annual  Variations.  Figure  2  shows  the  daily  cycle  of  tem¬ 
perature  at  three  different  heights  over  a  desert  area  for  summer  and  winter. 

The  difference  between  the  amplitude  of  the  January  and  July  temperature  curves 
show  the  specific  effect  of  solar  radiation  on  the  ground  temperature.  The  in¬ 
fluence  of  cloud  cover  on  the  daily  range  of  temperature  is  depicted  in  Figure 
3.  Both  figures  clearly  show  that  the  ground  temperature  responds  to  the  in¬ 
coming  and  outgoing  radiation  and  that  the  air  temperature  follows  the  ground 


Figure  2.  The  Dally  Cycle  of  Temper¬ 
ature  at  Three  Heights  (from  R. 
Geiger,  The  Climate  Near  the  Ground, 
Figure  377  harvard  University  tress, 
Cambridge,  Mass.  1950). 


Figure  3.  The  Dependence  of  the  Di¬ 
urnal  Temperature  Range  on  Cloudiness 
(from  R.  Qelger,  The  Climate  Near  the 
ground.  Figure  36,  Harvard  University 
frees,  Cambridge,  Mass.  1950). 
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temperature.  The  diurnal  range  of  temperature  is  smaller  as  the  distance  from 
the  ground  increases.  Needless  to  say*  these  patterns  are  familiar  to  the  prac¬ 
ticing  meteorologist  in  the  weather  station.  However,  it  is  important  to  be  not 
only  aware  of  the  general  behavior  of  the  temperature  at  a  given  level,  but  to 
understand  the  physical  basis  for  the  variation  and  its  implications  on  the  var¬ 
iation  of  the  temperature  gradient,  which  is  the  most  important  single  atmos¬ 
pheric  variable  affecting  the  dilution  of  air  pollution. 

Typical  vertical  temperature  profiles  of  a  land  station  of  the  mid-latitudes 
are  shown  for  specific  hours  in  Figure  4.  The  diurnal  variation  of  the  all- 
important  temperature  gradient  is  illustrated  by  Figure  5.  The  three  sections 
of  this  figure  also  include  the  effect  of  cloudiness  at  three  wind-speed  ranges. 
The  seasonal  variation  of  lapse  rates  at  a  southern  Idaho  station  is  presented 
in  Figure  6,  showing  the  temperature  difference  between  heights  Of  1.5  and  180 
meters.  Predictably,  the  steepest  lapse  rates  occur  in  the  middle  of  the  day  in 
summer.  The  maximum  average  inversions  also  occur  in  the  summer  months.  This 
might  seem  surprising,  but  it  is  a  direct  consequence  of  the  greater  diurnal 
range  of  ground  temperature  which  occurs  in  the  summer.  The  wintertime  inver¬ 
sions  are  more  persistent  than  those  of  the  summer,  frequently  lasting  for  days, 
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Figure  4.  Typical  Profiles  of  Temperature  in  the 
Lowest  1500  Feet  for  Various  Times  of  Day  (from 
Holland  [12] ) . 
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Figure  5-  Diurnal  Variation  of  Temperature  Difference 
Between  11  m  and  110  m  for  Different  Cloudiness  and 
Wind  Speed  [19] . 

tut  their  average  intensity  is  less.  Because  of  the  effect  of  the  lapse  rate 
on  eddy  conduction  of  heat,  it  is  the  prolonged  slight  inversions  associated 
with  strong  winds  in  winter  that  actually  result  in  the  greatest  downward  net 
transport  of  heat  in  the  atmosphere  near  the  ground. 

b.  The  Influence  of  the  Underlying  Surface.  So  far  we  have  discussed  the 
features  of  the  temperature  pattern  that  would  be  present  over  flat,  uniform  ter 
rain.  In  paragraph  2.b,  we  discussed  factors  that  would  influence  the  response 
of  the  ground  layer  to  the  radiation  balance.  Because  turbulent  mixing  in  the 
atmosphere  is  so  dependent  on  the  local  temperature  gradient,  we  must  now  ex¬ 
amine  the  effects  on  the  local  temperature  structure  produced  by  variations  in 
the  response  of  the  underlying  surface.  We  will  consider  specifically  the  ef¬ 
fects  of  soil  type,  including  enow  cover  as  a  special  case,  and  the  markedly 
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Figure  6.  Seasonal  Variation  in  the  Vertical  Temperature 
Difference  Between  1.5  n  and  120  m  (southern  Idaho)  [5]. 


Figure  7-  Diurnal  Range 
of  Temperature  In  Scan¬ 
dinavian  Summer  for  Two 
KindB  of  Ground. 


different  response  properties  of  bodies  of  water. 

The  most  important  differences  in  the  underly¬ 
ing  surface  are  its  density  and  heat  conductivity. 
The  way  these  properties  influence  the  ground  tem¬ 
perature  is  illustrated  in  Figure  7.  The  solid 
curves  show  the  diurnal  variation  of  temperature  in 
a  dry,  sandy  soil.  The  density  and  thermal  con¬ 
ductivity  of  this  type  of  soil  are  both  low  because 
the  "dead  air"  trapped  In  spaces  between  the  grains 
of  sand  reduces  both  of  these  properties.  The 
dotted  curves  show  the  temperature  range  expected 
for  a  densely-packed,  rocky  soil.  The  heat  pene¬ 
trates  farther  into  this  denser  type  of  soil  so 
that  the  range  of  temperature  at  the  surface  is 
smaller.  The  dampness  of  the  ground  surface  also 
influences  the  thermal  properties  considerably. 

If  water  is  trapped  within  the  spaces  instead  of 
"dead  air,"  conductivity,  density,  and  specific 
heat  are  all  increased. 


The  low  temperatures  experienced  at  night  over  a  snow  surface  are  well- 
known.  They  occur  because  the  snow,  especially  if  it  is  fresh,  has  a  much  lower 
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denolty  and  conductivity  than  any  type  of  bare  soil.  The  most  extreme  inversions 
In  the  lower  atmosphere  in  winter  are,  therefore,  found  over  snow  cover.  The 

ability  of  a  snow  cover  to  reject  the  heating  effect  of  the  sun's  rays  in  the 

daytime  ic  generally  overrated.  A  moment's  reflection  Will  bring  to  mind  the 
fact  that,  barring  unusually  strong  cold  adveotlon,  there  is  a  substantial  rise 
of  temperature  on  ciear  days  even  when  there  is  a  complete  oover  of  fresh  snow. 

This  rise  is  due  to  the  fact  that  some  10  or  15%  of  the  incoming  solar  radiation 

does  penetrate  into  the  top  few  inches  of  the  snow  oover  and  thereby  is  able  to 

warm  the  top  layers  noticeably  because  of  the  previously-mentioned  low  density 

*  ' 

and  conductivity  of  the  snow  cover.  The  effect  of  a  snow  cover  may  then  be 
summarized  as  follows.  At  night,  a  snow  aover  produces  extreme  inversions.  By 
day,  a  snow  cover  reduces,  but  does  not  completely  eliminate,  the  superadlabatlc 
lapse  rate  that  normally  occurs  near  the  ground.  '  i 

Snow  cover  represents  something  near  the  extreme  condition  of  low  conduc¬ 
tivity.  Open  bodies  of  water  represent  the  opposite  extreme.  The  one  property 
of  water  that  distinguishes  it  most  from  solid  surfaces  is  that  water,  being  a 
fluid,  can  transfer  heat  in  the  convective  mode,  which  is  usually  much  more  ef¬ 
ficient  than  heat  transfer  in  solids.  Any  body  of  water  that  is  more  than  a 
puddle  will  have  currents  set  up  in  it  by  the  stress  of  the  wind  ana  these  cur-  . 
rents  result  in  very  effective  mixing  so  that  any  temperature  changes'  occurring 
are  distributed  through  a  deep  layer.  Of  secondary  importance,  but  still  sig¬ 
nificant,  is  the  fact  that  incoming  radiation  from  the  sun  penetrates  into  a 
fairly  deep  layer  of  water  instead  of  all  being  absorbed  at  the  surface.  These 
two  differences  result  in  smaller  seasonal  and  much  smaller  diurnal  changes  of 
the  air  temperature  at  the  surface  over  water  than  over  land.  The  result  is 
that  over  water  bodies  there  is  very  little  diurnal  air  temperature  change  un¬ 
less  it  is  caused  by  artvoction  from  a  nearby  land  area.  Over  an  inland  body  of 
unfrozen  water,  the  temperature  gradients  are  the  opposite  of  the  patterns 
shown  in  Figure  6;  the  most  unstable  periods  are  at  night  during  the  winter  and 
the  most  stable,  by  day  in  summer. 

c.  The  Influence  of  Topography  on  Temperature.  There  are  two  main  types 
of  topographic  effects  on  air  temperature  and  lapse  rates.  The  first  is  due 
simply  to  differences  in  elevation  and  the  second  due  to  differences  in  the 
slope  of  the  ground  surface. 

Considering  the  second  type  first,  the  fact  that  the  ground  is  sloping 
changes  the  hour  of  the  day  at  which  the  ground  surface  is  most  nearly  perpen¬ 
dicular  to  the  sun's  rays  and,  thus,  the  time  of  maximum  temperature  and  maxi¬ 
mum  lapse  rate.  For  example,  the  west  side  of  a  valley  will  warm  and  also  cool 
earlier  than  the  east  side  since  the  sun  travels  from  east  to  west.  Similarly, 
in  the  mid-latitudes  a  valley  that  opens  toward  the  south  will  generally  have  a 
better-developed  thermal  wind  system  by  day  than  one  that  opens  toward  the 
north. 

The  diurnal  range  of  temperature  at  the  floor  of  a  valley  is  larger  than  on 
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its  aidping  aides.  This  is  not  particularly 
surprising  since  frequently  superadlabatlc 
lapse  rates  prevail  by  day  and  inversions  by 
night.  However,  the  situation  is  not  quite 
that  Simple.  The  fact  that  the  ground  is  not 
level  means  that,  at  a  fixed  level,  the  air 
closer  to  the  sides  of  the  valley  cools  more 
at  night  than  the  air  above  the  middle  of  the 
valley  floor.  This  establishes  a  pressure 
gradient  which  causes  the  air  to  flow  down  the 
slopes  toward  the  valley  floor.  This  air  move¬ 
ment  as  a  type  of  wind,  will  be  disoussed  lat¬ 
er;  however,  its  effect  on  the  temperature  is 
to  cause  the  whole  valley  to  fill  with  cold 
air.  Thus,,  the  temperature  inversion  at  the 
floor  of  the  valley  is  not  nearly  as  intense 
as  are  those  on  an  open  plain.  On  the  other 
hand,  the  inversion  near  the  top  of  the  valley 
is  steeper  than  it  would  be  if  this  drainage 
flow  was  not  present.  The  resulting  distri¬ 
bution  of  temperature  in  this  situation  is  shown  in  Figure  8,  Even  with  the 
valley  circulation  effect,  however,  the  changes  caused  by  the  proximity  to  the 
ground  surface  are  not  completely  masked.  The  nighttime  and  daytime  Isotherm 
patterns  are  very  similar  but,  by  day,  the  warmer  air  is  on  the  valley  floor. 

4.  Local  Wind  Structure; 

a.  The  Mean  Vertical  Variation  of  the  Horizontal  Wind  Speed.  It  is  readily 
noted  that  the  wind  speed  one  foot  or  less  above  the  ground  is  quite  a  bit  weak¬ 
er  than  the  wind  at  the  same  time  near  shoulder  height.  This  decrease  of  the 
wind  speed  near  the  ground  is,  of  course,  a  direct  consequence  of  frictional 
drag  of  the  atmosphere  against  the  earth's  surface.  Well  above  the  surface,  the 
steady-state  wind  speed  is  determined  by  the  balance  of  pressure  forces  and  the 
Coriolis  force.  The  level  above  which  the  friction  of  the  ground  no  longer  in¬ 
fluences  the  wind  speed  in  direction  is  usually  termed  the  "gradient  wind  level." 
At  heights  between  the  ground  and  the  gradient  wind  level,  the  vertical  distribu¬ 
tion  of  wind  speed  is  determined  by  the  action  of  molecular  viscosity  at  the 
ground  level  and  by  "eddy  viscosity"  at  higher  levels  in  a  manner  completely 
analogous  to  the  way  the  distribution  of  temperature  is  controlled  by  the  three 
processes  of  heat  transfer. 

The  diurnal  pattern  of  wind  speed  is  controlled  largely  by  the  effect  of 
temperature  gradient  on  the  Intensity  of  turbulent  mixing  and,  hence,  on  the 
"eddy  viscosity."  During  the  day,  buoyancy  effects  greatly  enhance  the  turbu¬ 
lent  mixing  with  height  so  the  gradient  of  the  wind  speed  is  rapidly  reduced  to 
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Figure  8.  Patterns  of  Tem¬ 
perature  on  a  Hillside  (from 
R.  Geiger,  The  Climate  Near 
the  Ground,  Figure  110.  har¬ 
vard  University  Press,  Cam¬ 
bridge,  Mass.  1950) 
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Figure  9.  Average  Wind-Speed 
Profiles  Constructed  from 
Measurements  Made  During  a 
1-Year  Period  at  the  Oak  Ridge 
National  Laboratory  [12]. 


a  small  value.  At  night  the  effect  of 
buoyancy  tends  to  oppose  the  Increase  of 
turbulent  mixing  with  height  so  that  the 
gradient  of  the  wind  speed  is  more  uniform. 
Figure  9  shows  typical  changes  of  the  wind 
speed  distribution  between  day  and  night  at' 
a  mid-latitude  station.  The  most  obvious 
difference  is  the  well-known  decrease  of 
wind  at  night  in  the  layers  close  to  the 
ground,  but  equally  pronounced  is  the  in¬ 
crease  of  wind  at  night  at  heights  from  a 
few  hundred  feet  to  about  3000  feet. 

Another  parameter  affecting  the  distri¬ 
bution  of  wind  speed  is  the  roughness  of 
the  ground  surface.  A  rough  surface  causes 
more  turbulence  near  the  ground  than  does  a 
smooth  surface  so  that  vertical  mixing  Is 
more  effective.  This  enhanced  mixing  re¬ 
sults  in  a  smaller  value  for  the  gradient 
of  the  wind  speed  near  the  ground.  This 
smaller  gradient,  when  combined  with  the 
no-wind  condition  at  the  ground  surface, 
results  in  a  lower  wind  speed  near  the 


ground  than  over  a  smoother  surface,  and  a  more  rapid  increase  of  wind  with 


height  in  the  layer  above  the  anemometer  level  over  the  rough  surface.  These 


effects  are  illustrated  in  Figure  10. 


Figure  10.  Examples  of  variation  of  Wind  with  Height  Over  Different 
Size  Roughness  Elements  (figures  are  percentages  of  gradient  wind)  [4], 
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b.  Variation  with  Haight  of  tha  Horizontal  Wind  Dlraotion.  Anothar  featura 
of  tha  winds  in  the  surface  layer  is  well-known  to  tha  meteorologist.  That  Is 
that  the  winds  in  tha  surface  layer  do  not  blow  parallel  to  the  isobars,  but 
cross  tha  isobars  at  an  angle  from  high  pressure  toward  low  pressure.  Farther 
from  the  ground,  as  the  "gradient  wind  level"  is  approached,  the  wind  dlraotion 
becomes  more  nearly  parallel  to  the  direction  of  the  geostrophlc  or  gradient 
wind.  The  importance  of  this  turning  of  the  wind  with  altitude  is  that  the  net 
drift  of  a  cloud  of  pollutants  will  be  in  a  direction  between  that  of  the  geo- 
atrophic  wind  and  that  of  the  actual  surface  wind,  provided  the  pollutants  are 
being  mixed  through  a  fairly  deep  layer. 

c.  The  Influence  of  Geographical  Features  on  the  Local  Wind.  There  are  two 
primary  ways  in  which  the  topography  can  modify  the  local  wind  regime.  The 
first  of  these  is  through  the  thermal  effects  of  different  surfaces.  The  second 
is  the  dynamic  effect  of  hills  or  mountain  ranges  on  the  general  flow  pattern. 

A  common  example  of  the  first  type  is  the  "sea  breeze."  Because  of  its  high 
conductivity  and  specific  heat,  a  body  of  water  changes  temperature  during  the 
day  much  less  than  nearby  land  areas.  This  difference  of  temperature  of  the 
surface  in  turn  causes  a  difference  in  the  density  of  the  air  over  the  surface 
such  that  the  pressure  over  water  is  higher  by  day  and  lower  by  night  than  that 
of  adjacent  land  areas.  This  pressure  difference  causes  the  wind  to  have  a  com¬ 
ponent  directed  from  the  water  toward  the  land  during  the  day  and  from  the  land 
toward  the  water  at  night. 

Another  example  of  a  thermally -Induced 
wind  system  develops  in  mountainous  terrain. 

As  discussed  previously,  when  the  atmosphere 
cools  at  night,  tne  air  next  to  the  sloping 
sides  of  the  valley  is  more  dense  than  the 
air  in  the  center  of  the  valley  at  the  same 
elevation,  this  causes  the  air  on  the  slopes 
to  drain  toward  the  valley  center.  This 
drainage  ultimately  results  in  a  general  flow 
toward  the  mouth  of  the  valley.  Gravity-flow 
systems  of  this  type  will  usually  extend  to 
the  height  of  the  adjacent  ridge  tops.  Fig¬ 
ure  11  depicts  the  wind  regime  that  results 
from  a  typical  combination  of  valley  wind  with  a  light  prevailing  wind.  The 
opposite  thermal  pattern  develops  during  the  daytime  and  there  is  a  tendency 
for  upslope  and  up-valley  winds.  However,  up-valley  winds  are  not  usually  as 
distinct  because  the  greater  turbulent  mixing  during  daytime  destroys  the  small 
local  density  gradients  that  maintain  these  local  circulations. 

We  turn  now  to  the  strictly  dynamic  influences  of  topography.  Perhaps  the 
most  dramatic  example  of  the  influence  of  topography  on  the  general  flow  is  the 
mountain-wave  phenomenon,  which  often  produces  severe  updrafts  and  downdrafts  on 

15 


Figure  11.  Combination  of  Val¬ 
ley  Wind  with  a  Light  Prevail¬ 
ing  Wind  [20], 
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Plgure  12.  Flow  Perpendicular  to  a  Set  of  Ridgelines. 

the  lee  side  of  mountain  ranges.  The  same  phenomenon*  in  miniature*  occurs  in 
the  vicinity  of  any  ridge -valley  system.  While  these  local  circulations  on  the 
smaller  scale  do  not  often  constitute  a  hazard  to  aircraft,  they  can  drastically 
alter  the  pattern  of  spreading  of  any  concentration  of  airborne  material.  It  is 
not  too  uncommon  for  the  local  wind  in  a  valley  to  be  completely  reversed  from 
the  prevailing  flow  (see  Figure  12).  However,  if  the  prevailing  winds  are  not 
exactly  perpendicular  to  an  elongated  valley*  a  similar  phenomenon  known  as 
channeling  results.  This  is  illustrated  by  Figure  13- 

Of  all  the  local  wind  circulations*  only 
the  sea  breeze  has  been  successfully  ana¬ 
lyzed  in  a  quantitative  fashion  that  permits 
accurate  predictions.  Yet  the  micrometeor¬ 
ologist  must  be  aware  of  the  qualitative 
changes  that  are  caused  by  the  other  topo¬ 
graphical  influences  so  he  can  visualize 
the  diffusion  and  transport  patterns  that 
result  from  topographical  modification  of 

the  broad-scale  wind  patterns.  Figure  13.  Channeling  of  Wind 

by  a  Valley  [20], 
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Chapter  3 

BASIC  CONCEPTS  IN  DIFFUSION 


1.  General. 

In  Chapter  1  we  mentioned  turbulent  mixing  as  a  means  of  diluting  the  con¬ 
centration  of  pollutants  or  toxic  matter  in  the  atmosphere.  This  turbulent 
mixing  was  called  eddy  conductivity  in  Chapter  2  when  explaining  the  distribu¬ 
tion  of  temperature  near  the  ground.  In  the  present  chapter  we  return  to  the 
process  of  diffusion  of  matter  in  the  atmosphere.  The  fundamental  theoretical 
problem  of  diffusion  is  to  relate  the  spread  of  a  cloud  of  material,  and  the 
consequent  reduction  in  concentration  of  the  suspended  matter,  to  the  turbulence 
of  the  wind.  The  practical  problem  is  to  find  means  of  estimating  the  rate  of 
diffusion  from  simultaneous  measurements  of  purely  meteorological  variables, 
such  as  wind,  temperature,  and  their  fluctuations.  In  this  chapter  we  also  pre¬ 
sent  some  fundamental  definitions  and  elaborate  somewhat  on  the  mechanism  of 
diffusion  so  that  the  relationship  between  diffusion  and  the  more  directly 
measurable  meteorological  parameters  becomes  clearer. 

2.  Basic  Concepts  and  Definitions. 

Some  of  the  basic  characteristics  of  the 
diffusion  problem  can  best  be  Illustrated  in 
the  context  of  the  effect  of  specific  types 
of  sources  of  contamination.  We  will  assume 
that  the  toxic  or  harmful  material  is  in  the 
form  of  gases  or  very  small  particles  that 
faithfully  follow  all  of  the  complex  motions 
of  the  atmospheric  whirls  and  eddies. 

a.  Instantaneous  Point  Source.  No 
source  is  completely  instantaneous  and  the 
point  source  concept  is  an  idealization,  but 
this  definition  is  a  useful  approximation 
to  the  release  of  harmful  toxic  matter  in 
a  short  "puff'1  or  an  explosion.  A  small 
spill  of  highly  volatile  material  would  also 
be  classed  as  an  instantaneous  point  source. 

One  can  easily  imagine  three  different  types 
of  efi^ects  of  atmospheric  turbulence  depend¬ 
ing  on  the  relative  sizes  of  the  initial 
cloud  and  the  eddies  in  the  atmosphere.  These  are  illustrated  in  Figure  14.  If 
the  motions  in  the  atmosphere  consisted  only  of  a  mean  wind  and  turbulence  on  a 
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Figure  14.  Idealized  Diffu¬ 
sion  of  a  "Puff"  [20]. 
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scale  much  smaller  than  the  size  of  the  "puff/*  the  situation  might  appear  as  in 
Figure  l4aj  the  cloud  drifting  with  the  wind,  but  being  enlarged  by  the  turbu¬ 
lence  (small  arrowi).  On  the  other  hand,  if  all  the  eddies  were  larger  than  the 
cloud,  the  situation  would  be  more  like  Figure  14b.  The  third  possibility, 
eddies  of  the  same  size  as  the  cloud,  would  produce  the  complex  result  of  Fig¬ 
ure  14c.  Actually,  turbulence  in  the  actual  atmosphere  consists  of  eddies  of 
all  different  sizes  so  that,  in  reality,  all  three  effects  are  combined.  For  an 
Instantaneous  joint  source  the  amount  of  material  Injected  into  the  atmosphere 
is  usually  given  as  the  total  quantity  of  material  released  (grams  or  milligrams), 
and  the  distribution  is  usually  given  in  terms  of  exposure,  which  is  the  inte¬ 
gral  of  the  concentration  as  the  cloud  passes  a  point. 

b.  continuous  Point  Source.  This  idealization  has  received  the  greatest 
amount  of  attention  in  the  literature.  Good  approximations  of  it  are  releases 
from  stacks  and  accidental  releases  caused  by  pipeline  cracks  or  spills  that  do 
not  quickly  evaporate.  The  quantity  of  material  released  by  a  continuous  point 
source  is  expressed  in  terms  of  rate  (e.g.,  kilograms  per  hour).  If  the  re¬ 
lease  is  prolonged  and  the  character  of  the  mean  wind  and  turbulence  do  not 
change,  it  is  appropriate  to  talk  about  the  average  concentration  of  the  con¬ 
taminating  material  at  a  given  point  downwind  from  the  Bource.  The  concentra¬ 
tion  will  fluctuate,  but  a  mean  concentration  over  a  chosen  sampling  period  is 
meaningful.  Figure  15  shows  how  the  irregular  fluctuations  smooth  out  with 
time  to  form  predictable  average  concentrations. 
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Figure  15.  Instantaneous  and  Average  Aspects  of  a 
Smoke  Plume  [20]. 


c.  Continuous  Horizontal  Line  Source .  This  problem  received  considerable 


attention  by  early  British  researchers  (e.g.,  Sutton)  in  connection  with  chemi¬ 
cal  warfare  developments.  This  type  of  model  is  now  finding  use  In  the  analy¬ 
sis  of  urban  air  pollution.  Any  consideration  of  this  source  type  by  AWS  mete¬ 
orologists  will  probably  be  limited  to  supporting  the  use  of  Incapacitating 
gases . 
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d.  Continuous  Area  Source.  This  source  configuration  Is  usually  modeled 
by  considering  the  area  aa  a  series  of  continuous  line  sources  oriented  perpen¬ 
dicular  to  the  prevailing  wind.  It  Is  of  considerable  interest  m  urban  air 
pollution.  The  main  application  for  AWS  meteorologists  would  be  for  a  detailed 
analysis  of  the  hazard  associated  with  large  spills  of  slowly-evaporating  li¬ 
quids, 

e.  Instantaneous  Horisontal  Line  Source.  This  possibility  corresponds  to 
an  aircraft  spray  source.  It  has  been  studied  by  the  Army  in  connection  with 
CB  weapons.  The  analysis  is  similar  to  the  instantaneous  point  source  excapt 
that  the  concentration  in  the  crosswind  direction  is  uniform. 

f.  Instantaneous  Area  Source.  This  case  is  also  of  interest  mainly  for 
Army  support.  The  analysis  consists  of  superimposing  the  effects  of  a  large 
number  of  parallel  line  souroes. 

g.  Instantaneous  Vertical  Line  Source.  This  case  corresponds  to  the  trail 
of  exhaust  fumes  from  a  rocket  launching  and  is  therefore  of  considerable  im¬ 
portance  to  AWS  staff  meteorologists  at  the  missile  ranges.  However,  the  asso¬ 
ciated  diffusion  analysis  is  theoretically  the  most  complex,  so  only  the  diffu¬ 
sion  of  the  ground-level  portion  of  the  trail  has  been  studied. 

3.  Theoretical  Analysis  of  Diffusion. 

There  is  no  single  physical  model  capable  of  explaining  all  the  features  of 
turbulent  diffusion.  There  are  two  basic  approaches,  each  of  which  has  merits 
not  possessed  by  the  other.  One  relates  the  transport  to  the  gradient  of  the 
concentration  by  an  exchange  coefficient.  For  example,  in  the  case  of  heat 
transport,  the  exchange  coefficient  is  the  same  thing  as  the  eddy  conductivity 
discussed  earlier.  The  second  approach  considers  the  statistics  of  the  air 
motions. 

a.  The  Gradient  or  "K11  Theory.  The  basic  hypothesis  is  that  the  rate  of 
diffusion  depends  in  a  linear  fashion  on  the  gradient  of  pollution.  This  ap¬ 
proach  was  first  proposed  by  the  German  physiologist  Adolph  Fick  to  explain  the 
diffusion  of  salt  in  a  liquid.  For  the  simplest  case  of  a  gradient  in  one  di¬ 
rection  and  uniform  concentrations  in  the  other  directions,  the  differential 
equation  describing  the  process  is 

Here,  q  represents  average  concentration,  t  is  elapsed  time,  and  x  is  distance 
in  one  direction.  K  is  the  coefficient  of  diffusivity.  The  solution  of  Equa¬ 
tion  (1)  is 

(2)  * =  jmz  exp 

Here  Q  represents  the  source  strength  (total  release  of  pollutant).  For  three- 
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dimensional  diffusion  from  a  point,  the  equivalent  equation  is 


Here,  x,  y,  and  s  are  the  distances  from  the  point  along  the  axes  of  a  conven¬ 
tional  right-handed  coordinate  system.  These  equations  are  straightforward  in 
application,  and  evaluation  of  air  pollution  would  be  comparatively  easy  if  a 
single  value  of  K  applied  to  the  whole  atmosphere.  As  we  shall  see  in  the  next 
two  subparagraphs,  the  actual  values  of  K  are  dependent  on  the  temperature  gra¬ 
dient;  therefore,  they  change  with  time  and  also  with  height.  Furthermore, 
they  also  depend  on  the  roughness  of  the  underlying  terrain  and  the  wind  speed. 
The  completely  general  differential  equation  would  be; 

<4)  af  "  lie  (*x  S)  +  ly  (*Sr  Sy)  +  *z  (Kz  H) 

Completely  general  solutions  to  this  equation  do  not  exist  in  analytical  (mathe¬ 
matical  function)  form.  Solutions  can  be  worked  out  on  electronic  computers 
for  particular  initial  concentrations  of  contaminant  and  particular  sets  of  as¬ 
sumptions  about  the  K  values. 

Even  if  the  above  equation  could  be  solved  in  a  convenient  manner,  there 
would  be  yet  another  objection.  The  K-theory  assumes  that  there  1b  no  transport 
in  a  given  direction  unless  the  concentration  has  a  gradient  in  that  direction. 
It  has  been  shown,  however,  (Smith  023])  that  if  there  is  a  vertical  wind  shear 
and  a  vertical  gradient  of  the  contaminant  concentration,  a  significant  hori¬ 
zontal  diffusion  of  the  contaminant  will  occur  even  if  there  is  no  horizontal 
gradient  of  the  concentration.  (This  is  quite  separate  from  the  horizontal 
transport  due  to  the  mean  wind.)  Thus,  the  complete  equation  of  diffusion 
would  have  to  contain  at  least  nine  additional  terms  to  describe  the  interac¬ 
tions  between  wind  shears  and  concentration  gradients. 

Having  expounded  on  the  difficulties  of  the  complete  K  theory,  it  is  only 
fair  to  say  that  many  simplifications  are  quite  plausible.  In  particular,  if 
the  vertical  dimensions  of  the  cloud  of  contaminant  are  small  compared  to  the 
distance  from  the  ground,  the  variations  of  the  K's  can  Bafely  be  neglected, 
leading  to  a  solution  of  the  form 

•i*  <KxVz)'*  [-W  (£  + .  f;  + 

When  such  simplifications  are  justified,  the  K  theory  gives  very  useful  results. 

b.  Statistical  Theory  of  Diffusion.  A  more  realistic  theory  of  atmospher¬ 
ic  diffusion  was  published  in  1921  by  the  noted 'British  fluid  dynamicist.  Sir 
Geoffrey  Taylor  Q331  •  The  most  useful  expression  of  this  theory  from  the  point 
of  view  of  understanding  turbulent  diffusion  has  been  given  by  Pasquill  [18]. 

It  relates  the  mean  square  distance,  y8,  traveled  by  a  particle  during  the  time 
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Interval,  t»,  to  the  power  spectrum  of  the  turbulent  velocity  fluctuation  in 
that  direction.  The  expression  is 

(5)  7“(ti)  -  l  ev(n)  pa.)8 

O' 

Here  v'8  is  the  variance  of  wind  speed  in  the  y  direction.  The  power  spectrum 
sv(n)  tells  how  much  of  the  turbulent  energy  v'8  can  be  regarded  as  equivalent 
to  an  oscillation  at  the  frequency  n.  The  term 

i^Sin  tt  nti  \* 

r^rsEirv 

has  a  maximum  value  of  unity  when  tj  ■  0  or  n  -  0.  Xt  can  be  seen  from  this 
that  the  value  of  the  Integral  in  Equation  (5)  depends  on  the  whole  speotrum  of 
turbulent  energy}  but,  for  the  larger  values  of  t*,  the  diffusion  is  dominated 
by  the  low-frequency  (slowly-varying)  components  of  the  turbulent  velocity. 

This  dominance  by  the  low-frequency  components  is  in  accord  with  our  intuition. 
If  the  motions  of  the  particles  are  mainly  at  low  frequencies,  the  particles 
will  undergo  large  excursions  before  changing  direction.  On  the  other  hand,  a 
particle  whose  motion  could  be  described  by  a  high-frequency  oscillation  might 
have  a  large  mean  kinetic  energy,  but  its  excursions  could  be  short.  Of  course, 
we  cannot  Identify  individual  particles  of  air  with  particular  frequencies  of 
the  power  spectrum}  nevertheless,  the  analogy  helps  us  visualise  the  informa¬ 
tion  contained  in  the  statistics. 

It  can  be  shown  that  this  statistical  theory  predicts  diffusion  equivalent 
to  a  "K"  that  increases  linearly  with  time  of  travel  while  close  to  the  source. 
Thus,  the  distribution  of  smoke  from  a  stack  or  other  point  source  is  described 
at  first  by  a  !,K"  that  increases  linearly  downwind.  For  values  of  travel  time 
longer  than  five  or  six  minutes,  the  statistical  theory  gives  almost  the  same 
result  as  K- theory. 

4.  The  Dependence  of  Diffusion  on  the  Temperature  gradient. 

Whether  the  diffusion  process  is  described  in  terms  of  an  exchange  coef¬ 
ficient,  a  mixing  length,  or  an  autocorrelation  function,  atmospheric  diffu¬ 
sion  is  enhanced  by  an  unstable  lapse  rate.  We  were  implying  K-theory  for 
heat  transfer  when  we  discussed  the  concept  of  eddy  conductivity  in  Chapter  2 
and  the  dependence  of  the  eddy  conductivity  on  lapse  rate.  The  reason  is 
very  simple.  When  the  lapse  rate  is  unstable,  a  particle  displaced  verti¬ 
cally  will  experience  a  gravitational  force  tending  to  move  it  further  in  the 
direction  of  the  displacement.  Since  the  force  is  exerted  along  a  direction 
in  which  the  particle  is  already  moving,  the  kinetic  energy  of  the  particle 
is  increased.  Since  motion  in  all  three  directions  is  linked  by  the  equation 
of  continuity,  the  increased  kinetic  energy  of  the  vertical  component  of  turbu¬ 
lence  results  in  an  Increase  of  the  horizontal  turbulence  as  well.  The  reverse 
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situation  prevails  when  the  lapse  rate  is  stable.  Energy  Is  removed  from  the 
vertical  component  because  the  buoyancy  (gravitational-density)  forces  are  op¬ 
posed  to  the  direction  of  motion.  Again,  the  reduction  of  energy  of  the  verti¬ 
cal  component  causes  a  reduction  of  energy  in  the  horizontal  components. 

The  high  degree  of  correlation  between  lapse  rate  and  the  intensity  of  tur¬ 
bulence  can  be  seen  by  comparing  Figure  16  with  Figure  6.  Figure  6  shows  the 
seasonal  and  diurnal  variation  of  lapse  rate,  while  Figure  16  shows  the  seasonal 
and  diurnal  variation  of  oA,  the  standard  deviation  of  the  horizontal  wind  di¬ 
rection,  an  indicator  of  the  intensity  of  the  lateral  component  of  turbulence. 


Figure  16.  Sea¬ 
sonal  and  Diur¬ 
nal  Variation  of 
the  Standard  Di- 
viation  of  the 
Horizontal  Wind 
Direction  at 
75  m.  (Esti¬ 
mated  from  Fig¬ 
ure  2.33  [20].) 
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A  slightly  more  elegant  ex¬ 
planation  of  the  effect  of  lapse 
rate  on  diffusion  is  provided 
by  the  statistical  theory.  As 
we  saw  in  paragraph  3b,  diffu¬ 
sion  depends  on  the  energy  in 
the  low-frequency  portion  of 
the  power  spectrum.  The  char¬ 
acteristic  time-scale  of  con¬ 
vective  motions  in  the  atmos¬ 
phere  is  from  one  to  several 
minutes.  Therefore,  convec¬ 
tive  motions  introduce  large 
amounts  of  energy  into  the  low- 
frequency  portions  of  the  spec¬ 
trum,  and,  thus,  the  integral 
in  Equation  (5)  will  be  large. 

In  contrast,  at  night  low- 
frequency  energy  is  selec¬ 
tively  absorbed.  Figure  17 
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Figure  17.  Composite  Spectra  of  Vertical 
Velocity  at  Brookhaven,  Long  Island  [16]. 
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chows  typical  turbulence  spectra  associated  with  various  lapse  rates.  The  curve 
labeled  "group  (a)n  is  for  quite  unstable  conditions.  The  curve  labeled  "group 
(o)"  is  for  near  adiabatic  lapse  rates. 

The  effeot  of  lapse  rate  on  diffusion  is  illustrated  by  some  typical  pat¬ 
terns  followed  by  the  effluent  from  a  smokestack,  which  have  now  become  classic 
in  literature  of  air  pollution.  These  patterns  and  their  relation  to  stability 
are  illustrated  in  Figure  18  and  discussed  below. 


Figure  18,  Effect  of  Lapse  Rate  on  Vertical  Diffusion  [20]. 

a.  Fanning.  This  condition  occurs  with  stable  conditions.  The  low- 
frequoncy  components  of  the  turbulence  are  suppressed,  especially  in  the  verti¬ 
cal,  so  the  vertical  spread  of  the  plume  is  very  small.  The  plume  may  show  a 
slow  meandering  of  direction.  With  respect  to  a  long-term  mean  wind,  this  me¬ 
andering  produces  an  effective  spreading  of  the  plume  which  is  much  greater 
laterally  than  vertically. 


b.  Fumigation.  There  is  a  brief  period  after  sunrise  when  the  lapse  rate 
is  stable  aloft  but  unstable  near  the  ground.  When  the  top  of  the  unstable 
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layer  reaches  the  height  of  the  stack,  the  diffusion  toward  the  ground  can  be 
large  while  the  upward  diffusion  remains  small.  A  high  concentration  of  pollut¬ 
ant  can  then  suddenly  be  brought  to  ground  level. 

c.  Looping  illustrates  the  muoh  larger  scale  of  tho  turbulent  eddies  under 
unstable  conditions.  Diffusion  is  very  rapid  so  that  points  away  from  the  axis 
of  the  plume,  but  near  the  source,  experience  a  greater  concentration  than  under 
more  stable  conditions.  On  the  other  hand,  further  downwind,  concentrations 
are  smaller  than  with  more  stable  lapse  rates. 

d.  Coning  is  characteristic  of  nearly  equal  intensities  of  turbulence  in 
the  vertical  and  horizontal  directions.  This  is  characteristic  of  temperature 
gradients  between  isothermal  and  adiabatic.  The  scale  of  the  eddies,  or  the 
low-frequenoy  energy  of  the  turbulence,  is  seen  to  be  less  then  experienced 
with  looping,  unstable  conditions,  but  more  than  in  the  stable  conditions  which 
produce  fanning.  Coning  is  characteristic  of  strong  winds  and  cloudy  skies. 

e.  Lofting  is  the  opposite  of  fumigation.  When  on  inversion  is  present  be¬ 
low  the  stack,  but  not  above,  the  preferred  direction  of  the  diffusion  is  upward. 
This  situation  is  very  desirable  from  the  point  of  view  of  avoiding  hazardous 
exposure  at  the  ground.  Of  course,  it  applies  only  to  elevated  sources  of  con¬ 
taminant  . 

5.  Variation  of  Diffusivlty  with  Height. 

Besides  being  dependent  on  lapse  rate,  the  diffusing  power  of  the  atmosphere 
also  varies  with  height.  Near  the  ground,  the  size  of  any  coherent  air  move¬ 
ment  will  be  limited  by  the  proximity  to  the  ground.  It  has  been  found  experi¬ 
mentally  that,  with  near-adiabatic  lapse  rates,  the  distribution  of  temperature 
and  wind  speed  are  explained  very  well  by  assuming  that  the  exchange  coefficient 
(K)  increases  linearly  with  height. 

As  would  be  expected,  the  statistical  theory  also  predicts  an  increase  of 
the  diffusivlty  with  height.  Investigations  of  the  spectra  of  turbulent  energy 
have  shown  that,  near  the  ground,  the  spectra  are  nearly  universal  functions  of 
the  non-dimensional  frequency,  f  ■  nz/U,  where  U  is  the  mean  wind  [L^l  •  The  for¬ 
mula  relating  spectra  under  a  change  of  frequency  variable  is  such  that 
sv(n)  ■  za /if  sv(f),  and  dn  «  U/z  df.  For  any  particular  height  and  atmospher¬ 
ic  stability,  the  turbulent  energy  is  proportional  to  the  mean  wind  speed,  e.g., 
v,TI/Ua  «  constant.  Therefore, 

*  a 

(6)  -Uztx8  [  sv(f)  (— g  yjfj"*)  df 

0 

Equation  (6)  shows  that,  if  the  turbulent  energy  is  constant  with  height,  the 
diffusivlty  will  increase  linearly  with  height  and  wind  speed.  Actually,  such 
a  quantity  as  v'B/Ua  decreases  with  height  in  the  manner  shown  by  Figure  19. 

When  the  atmosphere  is  stable,  the  decrease  of  turbulent  energy  with  height 
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Figure  19 •  Decrease  of  Lateral  Turbu¬ 
lence  with  Height  for  Different  Stabll 
ity  Conditions. 


is  very  marked.  In  addition,  un¬ 
der  stable  conditions ,  the  influ¬ 
ence  of  the  ground  surface  dis¬ 
appears  at  heights  of  only  a  few 
hundred  feet  and  the  speotra  no 
longer  obey  the  "universal"  rela¬ 
tionship  described  in  the  pre¬ 
ceding  paragraph.  Kalmal  and 
Haugen  [L4J  found  that,  in  stable 
conditions,  the  peak  of  the  ver¬ 
tical  velocity  speotium  becomes 
constant  with  height  rather  than 
being  inversely  proportional  tp 
height.  Thus,  in  stable  condi¬ 
tions,  the  maximum  vertical  dlf- 
fusivity  is  located  quite  dose 
to  the  ground.  The  decrease  of 
the  intensity  of  turbulence  with 
height  under  neutral  (adiabatic 
lapse  rates)  and  stable  condi¬ 
tions  eventually  dominates  the 
frequency  effeot  so  that  the  dlf- 
fusivlty  decreases  with  height  as 
the  gradient-wind  level  1b  ap¬ 
proached.  This  decrease  has  its 
counterpart  in  the  more  sophisti¬ 
cated  K-theories  which  have  been 
used  to  explain  the  wind  distri¬ 
bution  between  the  ground  and  the 
gradient-wind  level.  The  height 
dependence  of  K  from  one  of  these 
theories  [2]  is  shown  in  Figure 
20. 


Figure  20.  Typical  Variation  of  Verti¬ 
cal  Diffusivity  (K)  with  Height  (from  [2]). 

6.  The  Effect  of  Variation  in  the  Wind  Direction  with  Height. 

In  Chapter  2  we  mentioned  the  gradual  turning  of  the  wind  direction  between 
the  ground  and  the  gradient-wind  level.  When  a  diffusing  cloud  is  mixed  through 
a  deep  layer,  the  net  drift  of  the  cloud  will  be  in  a  direction  approximately 
halfway  between  the  surface  wind  direction  and  the  geostrophic  wind  direction. 
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Chapter  4 

ESTIMATING  DIFFUSION  AND  CONCENTRATION 
FROM  METEOROLOGICAL  OBSERVATIONS 


1.  General. 

In  this  chapter  we  discuss  the  practical  question  of  using  purely  meteoro¬ 
logical  observations  to  estimate  the  diffusion  and  calculate  relative  concen¬ 
trations  of  released  toxic  material.  Estimates  of  dispersion  of  toxic  matter 
can  be  made  by  two  distinct  methods,  air-tracer  surveys  end  calculation  from 
generalized  formulae. 

An  air-tracer  survey  consists  of  releasing  known  quantities  of  inert  gases 
or  small  particles  under  a  variety  of  weather  conditions  and  then  measuring  the 
quantity  of  tracer  material  that  arrives  at  various  locations.  The  observed 
concentrations  of  tracer  material  are  then  correlated  with  meteorological  pa¬ 
rameters  such  as  wind  speed,  temperature  lapse  rate,  wind  direction,  and  the 
variability  of  wind  direction.  From  these  correlations,  equations  are  developed 
to  express  concentration  directly  in  terms  of  the  meteorologies!  parameters. 
Examples  of  this  approach  are  diffusion  equations  developed  in  the  Mountain  Iron 
project  p.3]  and  the  Ocean  Breeze  and  Dry  Gulch  projects  [9J. 

Air-tracer  surveys  are  the  most  direct  and  accurate  method  of  estimating 
diffusion.  However,  this  approach  is  also  the  moat  expensive,  and  is  usually 
not  followed  until  a  definite  environmental  hazard  has  been  recognized  and  a 
practical  control  procedure  is  required.  Since  the  diffusion  equations  devel¬ 
oped  through  tracer  sampling  are  peculiar  to  the  location  where  the  survey  was 
conducted,  they  will  not  be  discussed  further.  The  balance  of  this  chapter  will 
present  more  general  methodB  which  can  be  used  to  make  preliminary  estimates  of 
diffusion  without  resorting  to  tracer  surveys. 

2.  Diffusion  Formulae  that  Assume  a  Gaussian  Distribution. 

Almost  all  of  the  formulae  commonly  used  to  describe  the  behavior  of  atmos¬ 
pheric  contaminants  assume  that  the  spatial  distribution  of  the  contaminant  oan 
be  described  by  the  bell-shaped  Gaussian  curve  used  in  statistics.  This  dis¬ 
tribution  is  appropriate  for  two  independent  reasons.  First,  a  distribution 
of  the  Gaussian  form  is  a  fundamental  solution  of  the  differential  equation  of 
diffusion.  Equations  (1)  or  (4),  if  the  K  or  K' s  are  constant.  Second,  a  well- 
known  law  in  statistics  statos  that  the  distribution  of  the  means  of  finite 
samples  of  random  variables  will  tend  to  the  Gaussian  form  as  the  number  of 
such  samples  tends  to  infinity.  This  principle  applies  to  concentration  meas¬ 
urements  of  air  pollution  which  are  averages  over  time  of  an  instantaneously- 
fluctuating  concentration. 
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a.  Sutton1 a  Dlf fusion  Modal.  One  of  tha  most  widely-used  Gaussian  models 
in  earlier  years  was  Sutton's  model.  While  it  was  first  proposed  as  an  appli¬ 
cation  of  Taylor's  statistical  theory  of  diffusion,  the  model  is  stated  in 
terms  of  parameters  that  resemble  the  diffusion  coefficients  of  the  K-theory. 

In  Chapter  3,  paragraph  3b,  we  introduced  the  ooncept  that,  after  diffusing  for 
a  time,  tx,  particles  would  have  a  mean-square  displacement  y*"  from  their  ini¬ 
tial  position.  Sutton's  theory  [22]  shows  that 

7  *  i  cy'  (uti)2"11 

where  Cy  is  a  "virtual  diffusion  coefficient,"  U  is  the  mean  wind  speed  and  n  is 
a  number  to  be  derived  from  the  power-law  expression  for  wind  profiles.  If  p  is 
the  power-law  exponent,  l.e., 

.  ,  1  ■ 

»<*>  -  u*  (!r)P 

then  n  *  2p/l+p.  For  the  commonly-used  1/7  power  law,  n  is  0.25.  For  an  in¬ 
stantaneous  point  source,  the  concentration  q  of  a  contaminant  is  given  by 


(7) 


q 


A. 


(TTUt)* 


c  e  c 

x  y  z 


exp 


[-  (ut) 


n-2 


This  is  identical  with  the  equation  at  the  end  of  Chapter  3,  paragraph  3a,  if 
n  »  1  and  Cx“  -  kK^U,  Cys  *  4Ky/u,  etc.  The  parameters  Cx,  Cy,  and  Cz  are  thus 
generalized  K-values.  In  the  formula,  t  is  measured  from  the  time  of  release 
and  n  is  a  parameter  associated  with  stability;  x,  y,  and  z  denote  distances 
from  the  center  of  the  cloud.  The  values  of  the  parameters  Cx,  Cy,  and  Cz  are 
to  be  determined  from  turbulence  statistics  by  the  formulae 


(I-nTCS-nT 


etc.  The  ratio  N/U  depends  on  height  above  the  ground  and  the  roughness  of  the 
ground  but  is  almost  independent  of  wind  speed.  N  is  a  parameter  Sutton  calls 
"macroviscosity."  The  parameter  U,!r  is  the  variance  of  the  x- component  of  air 
velocity.  As  mentioned  in  Chapter  3,  paragraph  5,  the  ratio  ir57ur  is  nearly 
constant  for  a  given  site,  height,  and  stability. 

The  steady-state  formula  for  a  continuously-emitting  point  source  at  ground 
level  is 


(8) 


q  (x,y,z) 


2Q 

TTC  C  0  2"n 
y z x 


In  this  equation  the  x-y-z  coordinate  system  remains  centered  on  the  source, 
rather  than  moving  with  the  mean  v/ind.  Two  other  parameters  derived  from 
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Sutton's  method  ere  the  maximum  concentration  at  the  ground  from  a  source  at 
height  h  above  the  ground. 


Uh  y 

and  the  value  of  x  where  the  maximum  ground  concentration  occurs. 


xmax 


The  Sutton  method  is  most  useful  for  nearly-adlabatic  lapse  rates  and  releases 
from  elevated  sources  where  the  vertical  variations  of  the  coefficients  are 
small. 

b.  More  Recent  Qausslan  Models.  While  Sutton' r  model  has  considerable 
theoretical  validity,  it  is  necessary  to  specify  four  parameters  n,  Cx,  Cy,  and 
Cjj,  in  order  to  use  the  formula,  all  of  which  are  dependent  on  height  and  sta¬ 
bility.  Actually,  a  Gaussian  distribution  is  completely  determined  by  its  mean 
and  its  standard  deviation.  The  modern  tendency,  therefore,  is  to  omit  the 
theoretical  considerations  and  to  simply  relate  the  standard  deviations  of  the 
contaminant  distribution  directly  to  the  meteorological  parameters.  Under  this 
concept,  the  form  of  the  distribution  for  a  continuous  source  located  at  ground 
level  is: 


<10>  5  -  whr  “P  [-  *  (£r  +  £•)! 

*  y  z 

In  this  equation,  as  in  Sutton's,  the  diffusion  along  the  mean  wind  axis  by  tur¬ 
bulence  has  been  neglected  in  comparison  with  the  transport  due  to  the  mean 
wind.  The  parameters  oy  and  oz  are  the  standard  deviations  of  the  Gaussian  dis¬ 
tributions  in  the  y  and  z  directions.  The  equivalent  equation  for  an  instan¬ 
taneous  point  source  at  ground  level  would  be: 


(11) 
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where  the  x,  y,  and  z  follow  the  cloud.  In  both  of  these  formulae  the  values  of 
°x*  °y»  and  az  are  not  constant,  but  increase  downwind  from  the  point  of  release. 
The  peak  concentration  at  the  center  of  the  cloud  decreases  downwind  while  the 
curve  representing  the  distribution  of  concentration  gets  broader.  On  the 
centerline  of  the  cloud,  the  values  of  y  and  z  are  zero;  therefore,  the  formula 
for  the  peak  concentration  Cp  becomes 


(12) 


CP“ 


The  maximum  ground-level  concentration  from  an  elevated  release  is 
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These  formulae  are  quite  praotical  when  the  values  of  ox,  o^.,  and  °z  are 
determined.  Methods  for  estimating  these  are  presented  in  the  following  para¬ 
graph. 

3*  Estimating  Diffusion  Parameters  from  Ordinary  Meteorological  Data. 

In  order  to  apply  equations  such  as  (10)  and  (12),  it  is  necessary  to  know 
0y  and  oz  as  a  function  of  distance  from  the  source.  Over  the  past  several 
years  there  has  been  a  continuing  effort  on  the  part  of  several  meteorological 
Agencies  to  determine  functional  relations  of  oy  and  os  to  distance.  In  spite 
of  this  continuing  effort,  there  is  no  agreement  on  a  single  "best  formula"  for 
the  growth  of  cy  and  ®z  downwind  as  a  function  of  normally- available  meteorolo¬ 
gical  information,  the  most  recent  studies  [25]  show  that  cA,  the  standard  de¬ 
viation  of  the  horizontal  wind  direction,  is  the  best  single  meteorological  pa¬ 
rameter  for  estimating  diffusion  from  continuous  sources.  This  parameter  is  a 
direct  measure  of  crosswind  turbulence  and  is  also  highly  correlated  with  the 
intensity  of  vertical  turbulence.  One  of  the  earliest  systems  for  relating 
diffusion  to  meteorological  measurements,  and  still  one  of  the  best  systems, 
was  developed  by  Cramer  [3]  based  on  measurements  of  erA  and  c^,  the  standard 
deviation  of  the  vertical  wind  angle.  These  parameters  are  difficult  to  calcu¬ 
late  without  the  use  of  analog  or  digital  computers,  so  other  systems  have  been 
developed  which  require  less  sophisticated  Information. 

The  system  used  by  meteorologists  supporting  the  Atomic  Energy  Commission 

and  the  National.  Air  Pollution  Control  Administration  was  originally  devised  by 

Pasquill  [17].  It  was  stated  originally  in  terms  of  the  height  and  angular 

spread  of  a  cloud  of  smoke.  The  relations  for  the  height  and  angular  spread 

were  redefined  in  terms  of  o  and  o  by  Gifford  [7]*  With  this  redefinition, 

y  z 

the  predictions  of  Pasquill  agree  with  Cramer  within  a  factor  of  two.  Pasquill* s 
graphs  (as  modified  by  Gifford)  relating  and  cz  to  distance  are  presented  as 
Figures  21  and  22.  Six  stability  categories,  A  through  P,  were  defined  by  Pas¬ 
quill  in  terms  of  wind  speed  and  a  qualitative  assessment  of  the  strength  of  the 
incoming  or  outgoing  radiation.  A  summary  of  rules  for  selecting  the  appro¬ 
priate  category  is  presented  in  Tables  1  and  2.  For  concentrations  on  the  plume 
axis.  Equation  (12)  indicates  that  Cp/Q  is  inversely  proportioned  to  Vay°z • 

Figure  23  depicts  CpU/Q  as  a  function  of  distance  for  the  six  categories. 

If  information  on  the  standard  deviation  of  horizontal  wind  direction  is 
available,  the  following  table  of  equivalence  (Table  3)  is  recommended.  The  val¬ 
ues  in  the  table  are  based  on  approximately  200  diffusion  trials  conducted  in 
the  Green  Glow,  Prairie  Grass,  Ocean  Breeze,  and  Dry  Gulch  projects  sponsored 
by  AFCRL,  and  projects  in  support  of  the  AEC  at  Idaho  Palls,  Idaho  and  Hanford, 
Washington.  As  with  Table  2, if  two  letter-categories  are  indicated,  one  should 
choose  a  value  halfway  between  the  two  curves. 
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A  rough  measure  of  oA  con  be  obtained  from  the  trace  of  wind  direction  from 
the  R02  wind  recorder  or  a  similar  instrument.  The  technique  consists  of  find¬ 
ing  the  extreme  range  of  wind  directions  during  a  30-minute  period  and  dividing 
by  eight  (15] •  The  method  is  very  crude  and  some  workers  prefer  to  divide  by 
only  five  or  six  instead  of  eight.  But  crude  as  it  Is,  the  method  is  more  ac¬ 
curate  than  assigning  the  Pasqulll  category  purely  on  the  basis  of  the  informa¬ 
tion  in  Tables  1  and  2. 

CAUTION:  Categories  A  and  B  should  never  be  used  for  nighttime  peek 
concentrations,  and  category  C  would  be  appropriate  at  night  only  when 
strong  gusty  winds  persist  following  a  cold  front  passage.  Values  of 
©A  at  night  larger  than  about  8°  over  flat  ground  or  14°  over  rough 
ground  indicate  fanning.  Short-term  peak  concentrations  will  then  be 
as  large  as  those  given  by  curve  F,  or  even  larger. 

Equation  (10)  can  be  evaluated  for  off-axis  points  using  Figures  23  and  24. 
For  the  Pasqulll  categories,  use  Figure  23  to  get  the  centerline  concentration 
appropriate  to  the  stability  category.  Figure  24  gives  the  value  of  the  expo¬ 
nential  term  for  lateral  distances  1  to  1000  meters  and  of  1  to  2000  meters. 
Figure  24  is  labeled  in  terms  of  general  numbers,  a  and  b,  because  it  can  also 
be  used  with  vertical  distances  and  <?z  to  get  concentrations  away  from  the 
ground.  Graphs  of  Uq/Q  for  all  points  on  the  ground-level  plane  have  been  pre¬ 
sented  in  [26]  based  also  on  the  Pasqulll  cry  and  ?z  values.  These  plots  are 
reproduced  as  Figures  25  through  30. 


Figure  24.  Graphical 
Aid  for  Computing 
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TABLE  1 

Net  Radiation  Claeses 


Elevation  of  Sun 

Cloud  Cover 

Night 

*  15* 

mu 

>  65° 

<  3/8 

3/8  or  4/8  at  any 

-2 

-1 

+i 

+2 

+3 

hat  or  broken  abv 
16000  ft  (inrl.  -©) 

©  abv  16000  ft  or 

-1 

0 

+2 

+3 

©  between  7000  ft 
and  16000  ft 

-1 

0 

0 

+1 

+1 

®  below  7000  ft 

0 

0 

0 

0 

+1 

©  below  7000  ft 

0 

_ _ 

0 

0 

0 

0 

TABLE  2 

Stability  Category  Based  On  Wind  Speed  and  Net  Radiation 


Surface  Wind  Speed 

Net 

Radiation  Class  (see  Table  1) 

Vsec 

knots 

+3 

+2 

+1 

0 

ma 

-2 

<  2 

<  4 

A 

A-B 

B 

D 

E 

2-3 

4-6 

A-B 

B 

C 

D 

E 

P 

3-5 

6-10 

B 

B-C 

C 

D 

D 

E 

5-6 

10-12 

C 

C-D 

D 

D 

D 

D 

>  6 

>  12 

C 

D 

D 

D 

D 

D 

TABLE  3 

Pasquill-Type  Classification  Bases  on 
Standard  Deviation  of  Wind  Direction 


— 

At  1/2  Mile  or 
Less,  Use  Type 

At  1.25  Miles, 
Use  Type 

At  5  Miles  or  More, 
Use  Type 

*  2.5° 

E-P 

E 

D-E 

5° 

D-E 

D 

C-D* 

H 

O 

o 

C-D* 

C* 

C* 

15° 

B* 

B-C* 

B-C* 

20° 

A-B* 

B* 

B-C* 

*  25° 

A* 

A-B* 

B* 

*  See  Caution  in  text. 
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Figure  25.  Isopleths  of 
WQ  on  the  Ground  for  a 
Ground-Level  Source, 

A  Stability. 


Figui 

qU/Q 


ire  26.  Isopleths  of 
...  on  the  Ground  for  a 
Ground-Level  Source, 

B  Stability. 


Figure  27.  Isopleths  of 
qu/Q  on  the  Ground  for  a 
Ground-Level  Source, 

C  Stability. 
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Pi gure  28.  Isopletha  of 
qU/Q  on  the  Ground  for  a 
Ground-Level  Source, 

D  Stability. 
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Figure  29.  Isopletha  of 
qlr/Q  on  the  Ground  for  a 
Ground-Level  Source, 

E  Stability. 


Figure  30.  Isopletha  of 
qU/Q  on  the  Ground  for  a 
Ground-Level  Source, 

F  Stability. 


4.  Non-Qauaslan  Equations. 

The  Gaussian  assumption  is  usually  retained  whenever  the  distribution  in 
more  than  one  direction  is  desired.  If  only  the  downwind  peak  concentration  is 
needed,  it  is  possible  to  drop  the  a  priori  functional  form  and  adopt  a  com¬ 
pletely  general  approach  in  relating  concentration  of  contaminants  to  meteoro¬ 
logical  parameters.  An  example  of  the  general  approach  is  the  technique  devel¬ 
oped  at  the  Air  Force  Cambridge  Research  Laboratories  in  connection  with  the 
Ocean  Breeze  and  Dry  Gulch  diffusion  programs  [9].  The  relevant  meteorological 
parameters  were  assumed  to  be:  the  mean  wind  speed  U,  the  downwind  distance  x, 
the  RMS  wind-direction  fluctuations  cA,  and  the  temperature  difference  AT  in  de¬ 
grees  Farenheit  between  six  feet  and  54  feet  (negative,  if  the  f4-ft  temperature 
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was  colder).  A  completely  general  equation  of  the  type 
Cp/Q  -  Kx^u'5  ce°  (AT  +  oonst)d 

was  fitted  to  the  Cp/Q  data  by  taking  logarithms  of  both  tides  and  using  multi¬ 
ple  regression  techniques  to  determine  best-fitting  values  of  K,  a*  b,  c,  and 
d.  Cp  is  the  peak  concentration  which  theoretically  corresponds  to  the  value 
of  q  on  the  centerline  of  the  Uaussian  model.  It  was  found  that  Cp/Q  was  inde¬ 
pendently  correlated  with  the  mean  wind  speed  U,  but  there  was  no  significant 
improvement  when  U  was  included  in  a  multiple  correlation.  The  resulting  equa¬ 
tion  was 

(14)  CpQ  -  0.00211x"1,96  oA“0,506  (AT  +  10)1**33 

AT  and  cA  are  strongly  correlated  with  each  other  as  indicated  in  Chapter  3, 
paragraph  4.  Therefore,  the  loss  of  accuracy  by  omitting  one  or  the  other 
would  be  small.  Since  0 A  is  hard  to  compute  accurately  without  digital  or  ana¬ 
log  computers,  a  simplified  equation  using  only  x  and  AT  was  developed  for  use 
at  the  Titan  II  launch  sites  [15] •  The  simplified  equation  is 


C p/Q  -  0.000175x"1-95  (AT  +  10)^’ 92 


U5) 

When  this  equation  was  tested  with 
a  set  of  independent  data  for  the 
same  locations,  the  predicted  con¬ 
centrations  agreed  with  the  ob¬ 
served  concentrations  65#  of  the 
time  within  a  factor  of  two  and 
within  a  factor  of  four,  94#  of 
the  time.  This  equation  is 
graphed  as  Figure  31. 

Most  applications  of  diffusion 
theory  to  safety  require  an  esti¬ 
mate  of  the  distance  beyond  which 
the  concentration  will  be  less 
than  a  specified  level.  For  this 
purpose.  Equation  (15)  can  be  con¬ 
verted  to  the  form 

(16) 

x  -  0.0388  (C  /Q)"0,513  (AT+10)2,53 

Mr 


ocmmm  mrmc*  <m) 
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Figure  31.  Peak  Concentration  vs  Down 
wind  Distance  from  Equation  (15). 
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In  deriving  these  equations,  data  from  three  different  sites  were  used,  one 
on  the  West  Coast,  one  on  the  East  Coast,  and  one  In  the  Great  Plains.  The  data 
from  the  three  different  sites  are  all  scattered  about  the  same  regression  lines} 
therefore,  the  relations  expressed  by  Equations  (14)  through  (16)  can  be  con¬ 
sidered  fairly  universal. 

5.  Instantaneous  Sources. 

The  Gaussian  equation  for  diffusion  from  an  instantaneous  source  was  given 
as  Equation  (11).  This  equation  requires  a  value  for  cx  as  a  function  of  down¬ 
wind  distance;  such  values  are  shown  in  Figure  32  for  three  broad  stability 
categories,  as  derived  by  the  US  Army  Munitions  Command  Operations  Research 
Group  [27].  The  values  of  are  larger  in  stable  conditions  because  true  down- 
wind  diffusion  depends  on  wind  shear,  which  is  largest  in  stable  conditions. 

The  best  values  of  cy  also  differ  from  their  values  for  continuous  sources.  The 
most  important  differences  are  that  oy  is  smaller  and  increases  less  rapidly 
with  distance  for  instantaneous  sources  than  for  continuous  sources.  Values  of 
oy  for  instantaneous  sources  are  given  in  Figure  33.  These  are  a  compromise 
between  results  of  tests  conducted  at  Dugway  Proving  Ground  [25]  and  values 
recommended  in  [20].  The  confidence  in  these  oy  values  is  less  than  for  contin¬ 
uous  sources. 


Technical  Report  214 


April  1971 


The  data  on  o  from  instantaneous 
cources  are  even  fewer*  Data  taken  from 
chemical -munitions  trials  at  Dugway  Prov¬ 
ing  Ground  185]  yield  estimates  for 
truly  ground*>level  releases  that  agree 
very  well  with  the  valuec  show  in  Figure 
22  for  continuous  soureos*  -  For  elevated 
sourest.  c*  values  given  in  [80]  suggest 
that  a  curve  between  C  and  D  of  Figure  22 
be  used  for  unstable  conditions,  a  curve 
between  D  and  E  for  neutral  conditions, 
and  a  curve  slightly  below  F  for  stable 
conditions .  Values  of  the  factor 
(Srr-o  d#}“*  for  instantaneous  , -round- 
level  releases  are  show**  in  Figure  34. 

6.  Compari son  of  Equations . 


So  far,  we  have  mentioned  two  dif¬ 
ferent.  systems  for  relating  air-pollution 
concentration?!  to  more  rCadily-available 
meteorological  data.  Neither  of  these  systems  can  be  considered  really  superior 
to  the  other.  Both  give  very  similar  answers  at  short  distances,  wher?  most  of 
the  data  have  been  gathered.  Unfortunately,  there  has  never  been  a  diffusion 
test  with  a  simultaneous  measurement  of  the  parameters  needed  to  tect  both  sys¬ 
tems.  In  most  cases,  the  choice  of  systems  is  dictated  by  the  type  of  data 
available .  e:i 


Figure  34.  (2t»  ayaz)~'m  tor  In¬ 

stantaneous  Relec.ee .  ' 


Actual  vertical  temperature  difference  measurements  are  seldom  available 
except  at  the,  Air  Force  test  ranges,.  Titan.  II  sites,  and  a  few  Army  installa¬ 
tions.  Nor  is  there  any  widely-accepted  system  for  predicting  these  vertical 
temperature  differences  from  other  meteorological  variables.  Thus,  for  most 
locations  the  only  available  tool  is  the  Paaoulll-Glfford  system  outlined  in 
paragraph  3  of  this  chapter  for  which  the  required  meteorological  data,  cloud 
cover,  wind  speed,  and  sun  angle.  Examples  of  calculations  with  this  system 
will  be  given  in  the  Appendix. 


7 .  Mixing  Depth;  The  Role  of  Elevated  Inversions. 

All  of  the  diffusion  equations  described  so  far  consider  only  the  effect  of 
the  lapse  rate  near  the  ground  in  determining  the  diffusion  (or  dilution)  from 
sources  at  or  near  the  ground. 

In  many  cases  the  lapse  rate  near  the  ground  will  be  isothermal  or  even 
adiabatic  with  a  pronounced  inversion  a  few  hundred  or  few  thousand  feet  above 
the  ground.  In  these  situations  the  diffusion  near  the  source  behaves  as  if 
there  were  no  inversion  present.  But  because  the  inversion  is  present,  the 
cloud  or  pollutant  cannot  grow  indefinitely  in  the  vertical  direction.  At 


-ywsr-'H  ^ 
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distances  far  enough  downwind  the  pollutant  is  uniformly  mixed  through  the 
layer  below  the  inversion  and  no  further  dilution  can  take  place  due  to  verti¬ 
cal  mixing.  The  cloud  does  continue  to  expand  horizontally  unless  the  lateral 
spread  is  also  confined  by  hills  parallel  to  the  flow. 

Some  .rules  of  thumb  for  the  effect  of  an  elevated  inversion  can  be  stated  as 
follows: 

a.  If  d  is  known  as  a  function  of  downwind  distance,  x,  determine  the  dis- 

Z 

iances  where  cr  «  H/2  end  where  0  *  2H,  H  being  the  height  of  the  inversion, 
z  z 

For  distances  such  that  cr2  <  H/2,  the  elevated  inversion  can  be  ignored. 

b.  For  distances  such  that  az  would  be  greater  than  2H,  assume  a  uniform 
distribution  in  the  vertical  and  use  the  formula 

(17)  q  -  — 3 - exp  (-y*/2o  *  ) 

VSt?  hu  oy  y 

c.  For  values  of  x  such  that  a  would  have  been  between  £h  and  2H,  the  con- 

z 

centration  can  be  estimated  by  drawing  a  straight  line  on  log-log  paper  between 
the  concentration  at  the  value  of  x  where  o_  would  be  &H  and  the  concentration 
at  the  distance  where  oz  would  be  2H. 

Unfortunately,  the  application  of  mixing-depth  relations  is  limited  to  the 
Gaussian  models  described  in  paragraphs  2  and  3  of  this  chapter.  The  statisti¬ 
cal  methods  of  paragraph  4  of  this  chapter  cannot  be  used  because  oz  is  not 
known  independently  as  a  function  of  distance.  The  inability  to  account  for 
mixing  depth  has  not  really  affected  safety  calculations  for  liquid-propellant 
spills  at  the  test  ranges.  The  mixing-depth  concept  is  important  for  air- 
pollution  problems  of  urban  areas  where  there  are  numerous  source's  of  pollution 
and  concentrations  at  the  "annoyance"  level  are  significant.  Spills  of  liquid 
propellants  of  the  magnitude  considered  in  test -range  safety  problems  would 
normally  be  diluted  to  safe  levels  in  less  distance  than  required  for  elevatpd  ‘ 
inversions  to  be  significant. 
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Chapter  5 

PROCESSES  OTHER  THAN  TURBULENT  MIXING  THAT  AFFECT  THE 
DISTRIBUTION  OF  TOXIC  MATERIAL 


1,  General . 

The  methods  presented  in  Chapter  4  assume  that  nothing  happens  to  the  ma¬ 
terial  other  than  dilution  by  turbulent  mixing.  This  is  not  entirely  the  case 
as  some  liquid  propellents  used  in  missiles  (also  some  CB  agents)  have  a  reac¬ 
tion  with  air  that  tends  to  diminish  their  concentration.  There  are  also  other 
processes  related  to  atmospheric  dynamics,  but  not  to  turbulence,  which  alter 
the  distribution  of  released  matter.  The  reaction  with  air  will  be  of  substan¬ 
tial  importance  for  some  missile  propellants  soon  to  be  tested.  This  type  of 
depletion  will  be  considered  first. 


2 •  Concentration  Changes  Resulting  from  a  Reaction  with  Air. 

Most  missile  liquid  propellants  of  the  hydrazine  family  react  with  air. 

This  reaction  has  been  recognized  for  quite  a  while  and  was  considered  when  de¬ 
termining  the  requirements  for  storage  facilities  but,  no  attention  has  been 
paid  to  the  effect  such  a  reaction  would  have  on  the  concentrations  in  the  at¬ 
mosphere  following  an  accidental  release.  In  some  cases,  as  with  unsymmetrical 
dimethyl  liydrazi.ne  (UDMH) ,  the  reaction  is  slow.  In  other  cases,  notably  mono- 
methylhydrazine  (MMH) ,  the  reaction  is  faster.  Within  half  an  hour  after  mixing 
with  adequate  amounts  of  air,  almost  half  of  the  MMH  iu  converted  to  hamless 
by-products,  which  substantially  reduces  the  downwind-hazard  distance.  The  re¬ 
action  may  be  even  faster  in  the  presence  of  certain  catalysts,  MMH  will  prob¬ 
ably  not  present  much  of  a  problem  for  the  AWS  staff  meteorologist  because  there 
are  no  current  USAF  plans  to  use  it  for  anything  except  small  rocket  motors  in 
deep-space  probes.  However,  the  reaction  of  MMH  provides  an  excellent  example  - 
to  illustrate  the  technique  of  calculating  modified  concentrations  for  reacting 
fuels . 

An  investigation  of  the  reaction  rate  cf  MMH  [29]  shows  that,  for  concen¬ 
trations  of  MMH  less  than  about  10#  by  volume,  the  rate  of  conversion  of  MMH  de¬ 
pends  on  the  amount  of  MMH  still  present.  That  is,  the  process  can  be  described 
by  the  first-order  differential  equation 


(18) 


dCMMH 

HT“ 


kCMMH 


in  which  is  the  concentration  of  MMH  and  k  is  the  proportionality  factor. 
The  solution  of  this  equation  is 
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in  which  CQ  Is  the  concentration  at  l00_ 

time,  t  »  0.  A  typical  plot  of  the  \ 

concentration  of  a  function  of  time  8  *°  \ 

is  shown  in  Figure  35 •  I  so ■  \ 

It  is  clear  that  this  behavior  5  \ 

S  40  ’  \ 

will  affect  the  pattern  of  diffusion  g  \ 

because  the  fundamental  concept  of  /  *°' 

the  statistical  theory  of  diffusion  — 

is  that  the  average  position  of  a  0  reaction  Tint 

particle  at  some  time,  t,  after  it 

Is  released  results  from  the  proper-  Figure  35.  Typioal  Concentration 

of  Rftmhiari  fu^iotiows  Ao-  of  Toxic  Substance  when  an  Initial'- 

ties  of  the  ensembied  functions.  Ac-  iy  Toxic  Sutstance  Decays  in  Air  to 

counting  for  the  details  of  the  dif-  Form  a  Harmless  Compound, 

fusion  j s  made  somewhat  more  compli¬ 
cated  by  having  to  account  for  the  changing  concentration  of  toxic  material  in 
each  "blob"  of  air,  but  the  most  Important  type  of  calculation,  the  peak  concen¬ 
tration  downwind,  can  be  corrected  very  simply.  As  a  first  approximation,  the 
process  of  dilution  due  to  the  chemical  reaction  can  be  considered  as  independ¬ 
ent  of  the  dilution  due  to  mixing.  In  that  case,  the  concentration  is  obtained 
by  multiplying  times  the  concentration  that  would  have  been  obtained  for 

a  nonreacting  substance.  Two  graphs  will  assist  in  making  this  calculation. 


Figure  35.  Typioal  Concentration 
of  Toxic  Substance  when  an  Inltial'- 
ly  Toxic  Substance  Decays  in  Air  to 
Form  a  Harmless  Compound. 


DECAY  RATE  (min-1)  DOWNWIND  DISTANCE,  X(nwtm) 

Figure  36.  Graph  for  Converting  De-  Figure  37.  Correction  Factor  for 

cay  Rate  and  Wind  Speed  to  K/U.  Concentration  of  Decaying  Pollutants. 

Figure  36  converts  values  of  the  decay  constant  k(min"1)  and  the  wind  speed  in 
knots  to  values  of  k/U  for  different  values  of  x  and  k/U.  Figure  37  gives  the 
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correction  factor,  e“kx//U,  for  different  values  of  x  and  k/U.  These  graphs  can 
be  used  to  correct  the  estimates  of  Cp/Cl  obtained  from  Pasquill's  curves,  or 
from  Equations  (14)  or  (15).  Paragraph  4,  Chapter  4  also  contains  an  equation. 
Equation  (16),  for  estimating  the  distance  beyond  which  the  peak  concentration 
would  be  less  than  a  prescribed  value.  The  corresponding  equation  for  a  BUb- 
stance  of  decaying  toxicity  would  be 

(19)  xp°*513hx/U  „  0. 0388 (Cp/q)"0, 513  (AT  +  10)2*53 

This  is  simply  Equation  (16)  with  x  replaced  by  xe0,313l5X//w.  The  value  of  x 
which  is  a  solution  of  Equation  (19)  is  naturally  smaller  than  the  value  ob¬ 
tained  from  Equation  (16).  Figure  38  furnishes  the  solution  for  Equation  (19) 
in  terms  of  the  solution  of  Equation  (16)  and  the  value  of  k/U  determined  from 
Figure  35. 


DISTANCE  CALCULATED  FROM  DIFFUSION  EQUATION  (m) 


Figure  38.  Revised  Hazard  Distances  for  Decaying 
Pollutants. 


So  far  we  have  discussed  the  case  of  a  toxic  chemical  that  reacts  with  air 
to  form  harmless  products.  The  reverse  possibility  exists,  i.e.,  a  /apor  that 
would  be  harmless  in  the  pure  state,  but  toxic  after  reacting  with  atmospheric 
oxygen  or  nitrogen.  The  equation  for  the  concentration  of  such  a  toxic  mate¬ 
rial  in  a  closed  chamber  would  be; 

(20)  C  -  Cf  (1  -  e“kt) 

where  C  is  the  concentration  of  the  new  toxic  compound  at  a  given  time  and  Cf  is 
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the  concentration  that  would  result  from  the  conversion  of  all  of  the  original 
chemical  to  the  new  compound.  A  typical  plot  of  concentration  versus  time  for 
such  a  reaction  is  shown  in  Figure  39. 


Figure  39*  Typical  Concentration 
of  Toxic  Substance  when  an  Initial¬ 
ly  Harmless  Substance  Reacts  with 
Air  to  Form  a  Toxic  Substance. 

Invoking  the  same  assumption  of 
independence  of  the  diffusion  and 
the  reaction  rate  which  we  used  pre¬ 
viously,  the  downwind  peak  concen¬ 
tration  would  behave  as  Indicated  in 
Figure  40.  This  behavior  is  quali¬ 
tatively  the  same  as  the  peak  con¬ 
centration  from  an  elevated  source j 
it  has  a  maximum  value  a  short  dis¬ 
tance  downwind  from  the  release 
point . 


Figure  40.  Generalized  Hot  of  Con¬ 
centration  Versus  Distance  for  a  Sub¬ 
stance  Having  Characteristics  Shown 
in  Figure  39- 


The  first-order  modification  to  the  downwind  concentration  can  be  found 
with  the  aid  of  Figure  37.  For  the  case  we  are  now  considering,  use  Figure  36 
to  obtain  k/U.  The  required  correction  factor  will  be  obtained  by  subtracting 
from  one  the  value  of  the  ordinate  from  Figure  37.  That  is,  if  the  value  from 
Figure  37  is  0.30,  the  required  correction  factor  is  0.70.  There  is  not  as 
much  point  in  calculating  a  correction  for  the  downwind  safe  distance  in  this 
case  because  the  large  distances  would  be  least  affected. 


Another  possible  concern  is  the  fact  that  some  materials,  such  as  hydro¬ 
chloric  acid  (HC1)  and  fluorine  compounds,  are  hygroscopic  and  could  be  diluted 
by  the  condensation  of  water  vapor.  There  is  some  disagreement  about  whether 
this  would  actually  reduce  the  toxicity.  In  such  cases,  the  Staff  Surgeon  or 
Staff  Bio-environmental  Engineer  should  be  responsible  for  setting  the  accept¬ 
able  concentrations  for  both  vapor  and  aerosol  forms.  If  the  aerosol  form  (as 
suspended  liquid  drops)  is  considered  to  be  less  harmful,  it  would  be  necessary 
to  estimate  the  rate  at  which  water  vapor  will  be  attracted  to  the  particles  of 
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toxic  cloud. 

3.  Other  Depletion  Processes: 

a.  Gravitational  Settling  and  Impacting.  The  exhaust  clouds  from  solid- 
propellant  rockets  contain  at  least  some  particles  that  are  large  enough  to 
settle  under  the  influence  of  gravity.  If  strong  acids  are  present  in  the  ex¬ 
haust,  smaller  particles  may  become  very  hygroscopic  and  accumulate  enough  water 
to  have  a  definite  fall  velocity.  In  general,  spherical  particles  larger  than 
20  microns  020]  will  have  a  tendency  to  settle  if  their  density  is  large  enough. 
For  liquid  water  drops  up  to  500  microns  in  diameter,  the  fall  velocity  of  the 
drops  is  aB  small  or  smaller  than  the  vertical  velocities  associated  with  day¬ 
time  turbulence.  The  appropriate  diffusion  model  for  this  case  is  a  combination 
of  normal  diffusion  and  settling.  A  procedure  for  revising  concentration  esti¬ 
mates  to  account  for  depletion  of  the  total  material  in  the  cloud  due  to  settl¬ 
ing  is  presented  in  Meteorology  and  Atomic  Energy,  1968  (2Cj|.  The  correction 
factors  corresponding  to  the  six  Pasquill-Oifford  stability  categories  are  shown 
as  functions  of  distance  in  Figure  41.  Fall  velocities  may  be  estimated  from 
Figure  42. 


Figure  4l.  Depletion  of  Pollu¬ 
tion  Concentration  Due  to  Gravi¬ 
tational  Settling  of  Particles 
(vf/U  «  .01  where  the  vf  is  the 

deposition  velocity). 


Figure  42.  Fall  Ve¬ 
locity  of  Spherical 
Water  Drops  as  a  Func¬ 
tion  of  Size. 


b.  Washout  by  Natural  Precipitation.  Washout  of  contaminants  from  the  air 
by  precipitation  is  of  considerable  importance  when  considering  the  concentra¬ 
tion  of  radioactive  material  in  the  atmosphere  because  it  represents  the  only 
way  other  than  gravitational  settling  by  which  radioactive  atoms  can  be  removed. 
For  that  reason  the  subject  of  precipitation  washout  is  Important.  A  useful 
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first  approximation  is  to  consider  that  the  amount  of  material  removed  from  the 
toxic  cloud  per  unit  time  is  proportional  to  the  concentration.  Then,  the  re¬ 
duction  in  concentration  downwind  ia  proportional  to  e_ltx//u  as  seen  in  the  pre¬ 
vious  section.  All  of  the  graphs  required  for  adjusting  the  downwind  concentra¬ 
tions  due  to  chemical  reaction  can  also  be  used  for  computing  washout.  The 
real  uncertainty  is  the  selection  of  a  proper  value  of  k  for  washout  (Section 
5.4  of  ($0).  The  washout  problem  is  relevant  to  two  types  of  situationn  in¬ 
volving  toxic  substances.  The  first  is  the  case  of  highly  soluble  gases  and 
the  other  is  when  the  toxic  material  is  in  the  form  of  solid  particles  too  small 
to  settle  out,  but  larger  than  one  micron  in  diameter.  It  is  not  expected  that 
staff  meteorologists  will  attempt  washout  calculations.  The  USAFETAC  should  be 
consulted  if  washout  calculations  are  needed  because  there  are  substantial 
problems  in  the  area  of  precipitation  climatology  involved  in  the  selection  of 
the  appropriate  k  values. 

4.  Characteristics  of  Real  Vapor  Releases. 

All  of  the  calculation  methods  presented  so  far  were  based  on  the  assump¬ 
tion  of  an  Instantaneous  point  source  or  a  continuous  point  source.  The  point- 
source  approximation  is  generally  acceptable  for  a  spill  of  liquids  covering  a 
fairly  large  area  as  long  as  the  concentration  calculations  are  not  needed  for 
downwind  distances  less  than  ten  times  the  diameter  of  the  spill.  Nevertheless, 
the  error  Involved  would  still  exceed  1$  for  a  distance  of  three  kilometers  if 
the  diameter  of  the  spill  area  were  50  feet.  The  error  in  the  point-source 
approximation  is  always  on  the  safe  side  for  locations  along  the  centerline  of 
the  cloud  if  the  calculation  is  made  assuming  that  all  of  the  toxic  vapor  orig¬ 
inates  at  a  point  in  the  center  of  the  spill  area. 

The  most  questionable  assumption  when  dealing  with  actual  vapor  releases  is 
the  assumption  that  the  source  strength  is  constant  for  15  minutes  or  more. 

This  assumption  is  especially  bad  in  the  case  of  either  compressed  gases  or 
cryogenic  fuels.  In  both  such  cases,  the  source  strength  decreases  rapidly  in 
the  first  few  minutes.  In  the  case  of  compressed  gases,  the  rate  of  release  is 
determined  by  the  internal  pressure  of  the  holding  tank.  The  tank  prese^fe  is 
reduced  with  time  not  only  because  the  gas  is  escaping,  but  also,  during  the 

I 

first  few  minutes  by  adiabatic  cooling  of  the  gas.  An  example.  Figure., 43  shows 

the  rate  of  escape  of  compressed  chlorine  through  a  1  1/3-inch  hole  in  a  55- 

* 

ton  tank  car  [13], 

The  rate  of  evaporation  of  a  spill  of  cryogenic  fuels  is  determined  primar¬ 
ily  by  the  temperature  of  the  ground  surface.  The  evaporation  cools  the  ground 
Just  as  the  temperature  of  a  wet-bulb  thermometer  is  reduced  by  evaporation. 
Eventually,  an  equilibrium  temperature  is  reached  at  which  the  heat  lost  from 
the  surface  due  to  evaporation  of  the  fuel  is  equal  to  the  upward  flow  of  heat 
by  conduction  through  the  soil.  A  recent  analysis  [30]  of  the  evaporation  of 
cryogenic  fuels  indicates  that,  during  the  cool-down  period,  the  evaporation 
rate  is  given  by 
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(21) 


«A 


Zi_i°  ihhe 

PL  AHl  v  nt 


where 


q  -  evaporation  rate 
A  per  unit  area 

T.  *  initial  tempera- 
1  ture  of  soil 

T  •  boiling  tempera- 
u  ture  of  cryogenic 

fuel 

k_  •  «•  thermal  conduc¬ 
tivity  of  soil 

C  ■  specific  heat  of 
8  soil 

Ps  -  density  of  soil 

PT  ■  density  of  fuel 
in  liquid  form 

Ahl  *  latent  heat  of 
evaporation  of 
fuel 


Equation  (21)  applies  from  short¬ 


ly  after  the  spill  to  the  time 
that  an  equilibrium  ground  tem¬ 


perature  is  reached.  The  change 
of  evaporation  rate  with  time  is 
illustrated  by  Figure  44.  The 


Figure  43.  Rate  of  Escape  of  Compressed 
Chlorine  Gtas  [13]. 


slope  of  the  line  on  a  log-log  plot  is  always  -1/2  because  the  denominator  of 


Equation  (21)  contains  the  square  root  of  time,  but  the  ordinate  value  depends 


on  the  physical  properties  of  the  soil  and  the  fuel  involved. 


The  non-steady  source  strengths  involved  with  actual  vapor  releases  result 


Figure  44.  Change  of  Evaporation  Rate  with  Time. 
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In  a  downwind  concentration  that  changes  with  time.  Typical  concentration  pro- 
files  at  two  different  times  are  shown  In  Figure  45  for  the  case  where  the  fuel 
of  Figure  44  is  spilled  Into  a  diked  enclosure ,  400  feet  on  a  side.  The  con¬ 
cave  shape  of  the  curve  reflects  the  decreasing  evaporation  rate.  For  any  par¬ 
ticular  time  after  the  spill,  there  would  be  two  downwind  concentration  maxima: 
one  close  to  the  spill  and  the  other  at  the  farthest  distance  downwind. 


1000  2000  3000  4000  8000  6000  7000  8000  9000  10,000 

DISTANCE  DOWNWIND  FROM  POOL  EDGE,  FT 


Figure  45.  Typical  Downwind  Profiles  of  Concentration  when  Evap¬ 
oration  Rate  Is  Not  Steady. 

For  fuels  which  have  a  boiling  point  near  or  above  normal  atmospheric  tem¬ 
peratures,  the  effect  of  evaporation  on  the  source  strength  is  much  less,  and 
it  is  possible  to  approximate  the  source  strength  by  a  steady  state.  There  is 
always  an  initial  wet-bulb  effect  but,  in  this  case,  the  change  of  evaporation 
rate  is  less  drastic  and  a  steady  state  is  reached  sooner.  The  prevailing  wind 
speed  determines  the  evaporation  rate  together  with  the  equilibrium  temperature 
The  most  appropriate  rate  of  evaporation  for  a  particular  fuel  should  be  ob¬ 
tained  from  operations  or  environmental  safety  personnel.  For  the  liquid  pro¬ 
pellants  used  in  the  Titan-series  rockets,  a  classified  report  by  Space  Tech¬ 
nology  Laboratories,  Inc.  [10]  contains  a  compilation  of  the  data  used  in  sped 
fying  the  operating-safety  distance  criteria  for  toxic  sources  at  the  Air  Force 
Eastern  Test  Range  and  the  Air  Force  Western  TeBt  Range.  That  report  is  con¬ 
sidered  the  best  currently-avallable  reference  on  the  relation  of  source 
strength  to  vind  speed,  propellant  temperature,  area  of  spill,  and  density  of 
the  liquid.  However,  since  it  Is  not  readily  available  to  most  AWS  personnel, 
we  strongly  recommend,  if  assistance  is  required  in  this  type  of  pollutant 
problem,  that  units  contact  USAFETAC  through  their  AWS  channels. 
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5.  Diffusion  from  Buoyant  Sources. 

So  far. we  have  considered  only  releases  at  ground  level  of  toxic  substances 
that  have  nearly  the  same  density  as  air.  This  limitation  was  appropriate  In 
safety  analysis  for  the  fuels  used  in  the  "core"  of  the  Titan-series  missiles. 
However,  there  are  now  rocket  motors  in  which  the  exhaust  products  are  toxic. 

The  classic  example  of  this  type  of  motor  is  the  solid-propellant  motor  with 
beryllium  additives. 

There  is  a  fundamental  difference  between  the  pollution  problem  of  missile 
liquid  propellants  and  the  pollution  problem  for  toxic  exhausts.  The  toxic  ex¬ 
hausts  are  very  hot  and,  therefore,  much  lighter  than  the  surrounding  air,  which 
causes  the  exhaust  to  rise  initially.  Opposing  this  tendency  for  the  exhaust  to 
rise  is  the  large  downward  velocity  of  the  gase3  in  a  normal  missile  launching. 
These  opposing  tendencies  result  in  the  formation  of  a  complex,  but  nearly 
spherical,  cloud  of  gases  at  the  launch  pad,  which  then  slowly  rises.  The  be¬ 
havior  of  this  cloud  for  Titan  III  boosters  has  been  analyzed  by  the  Aerospace 
Corporation  [8], 

a.  Traditional  Approaches.  The  toxic  hazards  associated  with  these  rising 
clouds  are  usually  evaluated  by  assuming  that  the  diffusion  pattern  at  any  dis¬ 
tance  downwind  is  equivalent  to  the  pattern  fr r,*\  a  non-buoyant  elevated  source. 
Near  the  source  this  equivalent  release  height  would  necessarily  become  higher 
and  higher  with  increasing  downwind  distance  as  the  exhaust  cloud  continues  to 
rise.  Eventually,  when  a  stable  layer  is  reached,  the  cloud-rise  will  stop. 

The  equivalent  release  height  would  then  b.'  constant  for  calculations  further 
downwind.  Fortunately,  for  most  toxic  haiard  calculations  the  equilibrium 
cloud-rise  is  achieved  within  a  few  minutes  and,  therefore,  a  constant  "effec¬ 
tive  height"  can  be  used.  Only  for  a  detailed  analysis  of  hazards  near  the  pad 
would  it  be  necessary  to  consider  the  change  of  "effective  height"  downwind. 

The  existing  data  on  the  rise  of  instantaneous  clouds  consist  mostly  of  ob¬ 
servation  of  "fireball"  clouds  from  small  nuclear-weapon  experiments  conducted 
by  the  Atomic  Energy  Commission.  Figure  46  shows  the  height  rise  as  a  function 


Figure  46.  Height  Rise  of  Instantaneous  Hot  Releases  [20]. 
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of  Qj/AS  where  is  the  heat  energy  released  and  66,  the  difference  between  the 
potential  temperature  at  the  top  of  any  ground-based  inversion  and  the  potential 
temperature  at  the  equilibrium  level.  The  solid  line  represents  the  equation 

Q  . ^ 

(22)  h  -  2.66  (•Q—p-lr/Sz) 

which  was  derived  from  small-scale  laboratory  experiments.  The  height  given  by 
Equatlor  (22)  may  be  combined  with  Equation  (9)  or  (13)  to  estimate  the  maximum 
concentration  at  ground  level.  However,  the  estimates  of  h  from  Equation  (22) 
are  not  too  reliable  and,  therefore,  the  concentration  estimates  would  be  only 
a  first  approximation. 

If  the  values  of  heat  released,  Qi,  are  at  all  large,  it  is  doubtful  whether 
the  maximum  ground-level  concentration  would  be  obtained  before  the  stability 
regime  changes  due  to  diurnal  effects.  The  real  ground-level  maximum  would  oc¬ 
cur  in  connection  with  a  process  known  as  "fumigation,"  which  was  mentioned 
briefly  In  Chapter  3. 

In  a  typical  fumigation  situation  the  toxic  products  are  initially  carried 
aloft  by  buoyancy  into  an  elevated  Inversion,  where  the  products  remain  sus¬ 
pended  with  little  diffusion.  After  sunrise  an  adiabatic  lapse  rate  1b  esti¬ 
mated  from  the  ground  up  to  the  base  of  the  stable  layer.  When  the  top  of  the 
adiabatic  layer  reaches  the  height  where  the  contaminants  are  located,  the  con¬ 
taminants  are  rapidly  mixed  throughout  the  adiabatic  layer.  The  concentration 
would  then  be  uniform  through  the  depth  of  the  remaining  inversion. 

Unfortunately,  there  are  no  tested  formulae  for  predicting  the  concentration 
due  to  buoyant  rise  and  subsequent  fumigation.  Two  formulae  have  been  proposed. 
One,  by  Briggs  [20]  for  a  steady  source  is 

(2os  r  „  _ _ 0^53 _ 

K  ^  PVb  a-*  (h  +  5.1  F*  s-*) 

where  q‘max  **  the  average  concentration  at  the  point  of  maximum  average  con¬ 

centration  (g/ma) 

Q.  *  rate  of  release  of  toxic  material 
F  ■  buoyancy  flux  given  by  P  ■  g(6T/Tg)  Vgr8 
g  -  gravitational  constant 

hfi  ■  height  of  stack 

r  *  stack  radius 

s  ■  stability  parameter  given  by  g/T  B0/8z 

T_  »  temperature  of  stack  gases 

S 

AT  -  temperature  difference  between  stack  and  environment 
V-  ■  vertical  velocity  of  stack  gases 
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The  other  formula,  developed  at  USAFETAC,  is  applicable  to  the  fumigation  from 
large  motors,  where  the  hot  exhaust-cloud  is  expected  to  be  distributed  through 
a  fairly  deep  layer  of  any  inversion.  The  USAFETAC  formula  i'. 


(84) 

where 


-  _  SAQ/At _ 

q  "  m  U(hc  +  h“/8) 

q  is  defined  as  in  Equation  (83) 

AQ  •  the  amount  of  material  that  stops  rising  between  any  pair  of  isen- 
tropic  surface^  0X  and  0*,  within  the  inversion 

At  *  the  time  Interval  over  which  the  temperature  at  the  ground  rises 
from  6X  to  0«  (minutes) 

h  «  the  height  above  ground  of  the  remaining  Jnve rsion  (meters) 

U  -  mean  wind  speed  at  the  inversion  (meters/minutes) 

oy  -  lateral  standard  deviation  of  cloud  at  the  point  and  time  where 
calculation  is  made 


The  summation  on  AQ  refers  to  the  amounts  of  material  previously  "fumigated" 
from  lower  levels  within  the  inversion.  The  equation  can  be  altered  to  fit  the 
case  near  a  shoreline  where  the  temperature  in  summer  increases  Inland.  In  that 
case,  UAT  should  be  replaced  by  Ax,  a  distance  over  which  the  temperature  rises 
by  the  required  amount. 

The  formula  requires  a  value  for  ey,  the  lateral  standard  deviation  of  the 
cloud.  An  estimate  of  cry  for  coastal  areas  which  fits  data  on  mesoscale  diffu¬ 
sion  given  by  [1]  and  also  those  data  presented  in  [20]  is  cr  ■  25  t,  where  t  is 
the  time  in  minutes  since  the  release  and  is  in  meters.  Also  fitting  the  data 
from  the  references  listed  above,  the  relation  Oy  »  JO  x  can  be  used,  where  x  is 
the  distance  inland  in  kilometers,  for  those  coastal  cases  where  the  temperature 
change  with  distance  is  more  significant  than  the  local  change  with  time.  These 
relations  probably  overestimate  the  rate  of  cloud  growth  with  time  or  distance. 
The  most  recent  experimental  work  gives  a  rate  of  spreading  that  is  proportional 
to  x0*®  or  t°  .  This  is  very  near  unity  and  substantial  differences  would  be 
noticed  only  for  very  large  travel  times.  Equation  (23)  contains  a  measure  of 
the  lateral  spreading  in  the  numerical  coefficient;  therefore,  no  further  esti¬ 
mates  of  spreading  are  needed  to  apply  that  formula. 


b .  Comment  on  Statistical  Techniques.  Conpletely  statistical  approaches 
to  the  calculation  of  downwind  hazard  distances  have  been  very  successful  for 
cold  ground-level  sources.  These  equations  were  discussed  in  Chapter  4.  The 
success  of  these  equations  has  made  it  very  tempting  to  try  the  same  techniques 
for  hot  sources.  However,  there  are  several  complications  inherent  in  trying 
to  derive  a  purely  statistical  one-step  prediction  model  for  buoyant  sources. 


These  difficulties  all  stem  from  the  fact  that  a  buoyant  source  is  effec¬ 
tively  an  elevated  source,  except  in  the  immediate  vicinity  of  the  launch  pad. 
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Prom  the  dieoussicn  in  the  preceding  paragraph,  we  know  that  the  ground-level 
concentration  will  have  two  maxima  and  two  minima.  Thn  primary  maximum . i?.  on 
the  launch  pad  itself.  The  concentration  will  decrease  rapidly  away  from  the 
launch  pad  because  the  exhaust  cloud  is  initially  rising  much  faster  than  it  la 
diffusing.  A  short  distance  downwind  the  cloud  center  will  have  leveled  off, 
and  from  that  point  onward,  the  ground-level  concentration  increases  until  a 
secondary  riiaxi’iium  is  reached  a  considerable  distance  downwind.  Finally,  beyond 
the  secondary  maximum  the  concentration  decreases  to, taro- as  the  distance  be¬ 
comes  very  large.  From  these  considerations  it  ,1s  evident  that  an  expression 
of  the  form  0p/Q  proportional  to  x"a  (a  continuously  decreasing  function  of  x) 
could  not  poscibly  describe  the  concentrations  for  all  values  of  downwind  dis¬ 
tance.  Furthermore,  the  dependence  of  Cp/Q  on  AT  or  o^  reverses  abruptly  at  the 
first  downwind  minimum  and  again  at  the  secondary  maximdfo.  In  the  zones  Where 
doncent ration  decreases  with  distance,  the  dependence  of  Cp/4  on  ST  end  would 
bs  qualitatively  similar  to  the  cold- source  situation.  The  dependence  of  Cp/Qi  ' 
on  these  parameters  is  reversed  in  the  region  between  the  first  downwind  mini¬ 
mum  and  the  secondary  snaximum.  Finally,  the  distances  to  the  first  minimum  and 
the  secondary  maximum  are  dependent  on  Qh  (the  heat  released)  and  AT. 

*■  v  ’  !  "  '  ■  '  ■  \  (' 

in  very  unstable  conditions  the  concentration  at  the  secondary  maximum  may 
not  be  much  larger  than  at  the  first  minimum.  A  scatter  diagram  of  concentra¬ 
tion  versus  distance  would  show  only  a  general  decrease  bf  concentration  with 
distance  and  no  particular  correlation  with  AT  or  erA.  This  in  exactly  the  re¬ 
sult  obtained  in  she  Sand  Storm  diffusion  program  [24],  which  was  carried  out 
only  wit a  unstable  conditions  at  Edwards  Air  Force  Base. 

.  '  .  v  i  1  ‘ '  ‘  ‘  . 

The  conclusion  that  prediction  equations  of  the  form  of  Equations  (14)  or 
(15)  will  not  work  for  all  possible  distances  does  not  mean  that  the  develop¬ 
ment  of  statistical  ohe-step  diffusion  equations  Is  impossible,  only  that  con¬ 
siderably  more  physical  insight  is  needed.  In  particular,  it  may  be  possible 
to  develop  one-step  Aquations;  for  limited  portions  of  the  downwind  hazard  area. 
Within  any  zone  where  the  concentration  is  increasing  or  decreasing  continuous¬ 
ly  with  distance,  a  simple  equation  of  the  type  used  in  Equation  (14)  would  work. 
If  it  is  possible  to  identify  certain  ranges  of  downwind  distance  for  which  the 
concentration  is  always  in  the  increasing  zone  or  always  in  the  decreasing  zone, 
separate  equations  could  be  developed  for  each  zone. 

The  only  parameter  needed  in  addition  to  those  in  Equation  (14)  is  the 
height  of  buoyant  cloud-rise.  This  height,  in  turn,  depends  on  and  AT.  To 
include  the  buoyant  rise  it  is  necessary  only  to  put  in  an  additional  factor 
Qhe  and  recalculate  the  multiple  regression.  The  bast-fitting  exponent  for  AT 
will  adjust  itself  to  reflect  the  role  that  AT  plays  in  determining  both  the 
cloud-rise  and  the  vertical  diffusion. 

c .  Continuous  Buoyant  Sources.  The  same  general  principles  apply  to  air 
pollution  emitted  from  smokestacks  or  longer-burning  fires.  The  initial  heat 
of  the  effluent  causes  the  plume  to  rise  and  become  equivalent  to  a  3ource 
higher  than  the  actual  height  of  emission.  The  traditional  approach  has  been  to 
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estimate  the  final  height  of  tho  plume  and  use  that  height  In  the  formula  for 
diffusion  from  an  elevated  source.  This  procedure  Is  acceptable  only  for  dis¬ 
tance  downwind  such  that  the  plume  has  nearly  reached  its  final  height. 

Until  recently,  this  restriction  on  downwind  distance  was  Ignored  In  most 
air-pollution  calculations.  However,  0.  A.  Briggs,  a  meteorologist  at  the  Oak 
Ridge  National  Laboratories,  has  given  a  simple  formula  for  "transitional"  plume 
heights  220] .  This  formula  for  the  change  of  effective  source  height,  Ah,  with 


distance  is 

(25)  Ah  -  2.0  F***  x*3  U"1 

where  U  »  mean  wind  at  Btack  height 

x  *  downwind  distance 

The  following  were  defined  in 
connection  with  Equation  (23) 

F  -»  buoyancy  flux, 

F  •  g ( AT/T )  Vgra 

g  -  gravitational  constant 

r  »  inside  stack  radius 

T_  *  stack  temperature 

AT  •»  temperature  difference 
between  stack  and  en¬ 
vironment 

Vg  »  stack  exit  velocity 

Once  the  buoyancy  flux,  F,  has 
been  calculated,  the  calculation  of  Ah 
is  easily  done  with  graphical  aids  as 
follows : 

(1)  Calculate  buoyancy  flux. 

(2)  Estimate  wind  speed  at 
release  level. 

(3)  Enter  Figure  47  with 
buoyancy  flux  as  the  abscissa  and  go 
up  to  the  intersection  with  the  wind 
value.  Then,  go  across  to  the  ordi¬ 
nate  to  find  the  value  of  2F^a  U""1 . 

(4)  Enter  Figure  48  with  the 
downwind  distance  as  the  abscissa  end 
go  up  to  the  sloping  line  correspond¬ 
ing  to  2F^  U”1 .  Frcn.  this  point,  go 


Figure  47.  First  Step  in  Evaluating 
Transitional  Plume  Rise. 


Figure  48.  Second  Graph  for  Trans i 
tional  Plume  Rise. 
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r cross  to  the  ordinate  to  get  the  value  Ah. 

(5)  Check  to  ree  that  the  final  plume  height  has  not  been  reached  (see 
following  paragraph).  ■' ,t  ,  ' 

It  Is  observed  that  all  plumes  and  buoyant  clouds  of  air  pollution  eventual¬ 
ly  lost  their  buoyancy.  The  prediction  ot  the  final  height  ia  especially  com¬ 
plicated  for  smokestacks  because  the  gases  start  out  with  some  verticil  momen¬ 
tum.  Consequently ,  a  large  number  of  different  formulae  have  appeared.  One > of 
the  most  successful  to  datye  was  developed  by  Holland  I.12J.  It  heads 

(26)  Ah  -  (i.S.  +  2.69  x  lO"8  p  -ft-qr-A  d) 

where  Ah,  u,  TR,Tg,  and  Vfi  are  defined  as  before,  ahu:  ,,  ''' 

d  w  Irveide  stack  diameter  r. 

p  -  atmospheric  pressure 

This  formula  is' now  found  to  underestimate  the  final  height.  Better  results 
are  obtained  by  multiplying  Ch  by, three. 
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Chapter  6 

THE  ROLE  OP  THE  STAFF  METEOROLOQIST  IN  AIR  POLLUTION  CONTROL 


1.  general. 

There  are  two  essential  steps  in  air  pollution  oontrol.  They  are:  (1)  re¬ 
ducing  emission  of  pollutants  where  practical  by  improved  engineering  and  (2) 
curtailing  operations  when  atmospheric  dispersion  is  not  active  enough  to  pro¬ 
duce  acceptable  pollution  levels.  A  corollary  activity  is  establishing  hasard 
areas  in  case  toxic  pollutants  are  involved.  Meteorological  support  is  needod 
for  operational  decisions  on  curtailing  operations  or  setting  up  potential 
hazard  areas,  and  for  planning  purposes  in  all  aspects  of  air  pollution. 

2.  Meteorological  Support  for  Planning. 

Meteorological  support  for  planning  purposes  is  just  another  kind  of  clima¬ 
tology  —  a  diffusion  climatology.  Normally,  new  facilities  will  be  sited  and 
constructed  so  as  to  minimize  the  possible  air  pollution  effects  of  any  planned 
or  accidental  release.  If  it  were  possible  to  find  a  perfect  site,  one  at 
which  pollution  concentrations  would  never  exceed  safe  levels  in  populated  areas, 
there  would  be  no  need  for  operationally-oriented  meteorological  support. 

The  staff  meteorologist  normally  reviews  plans  with  the  goal  of  determining 
what  sort  of  meteorological  support  will  be  required.  The  factor  of  possible 
air  pollution  should  now  be  included  in  this  review.  It  is  impossible  to  em¬ 
phasize  too  strongly  the  importance  of  examining  ell  aspects  of  Air  Force  re¬ 
search  and  development  programs  to  determine  whether  a  serious  environmental 
hazard  may  exist  routinely  or  accidentally.  While  the  recognition  of  environ¬ 
mental  hazards  is  normally  the  responsibility  of  bio-environmental  engineers, 
it  is  often  the  meteorologist  who  has  the  better  "feel"  for  those  combinations 
of  emission  and  meteorological  conditions  that  will  produce  dangerous  concen¬ 
trations  of  pollutants.  The  price  of  falling  to  detect  environmental  hazards 
until  late  in  the  development  program  is  often  a  costly  re-programming  of  the 
system.  If  the  environmental  hazard  is  recognized  early  in  the  development, 
meteorological-data  networks  and  control  procedures  can  be  developed  concur¬ 
rently. 

Mitch  of  what  has  just  been  said  pertains  to  such  pollutants  as  toxic  missile 
fuels  and  combustion  products.  But  there  is  also  a  need  for  meteorological  sup¬ 
port  In  planning  the  more  routine  air-base  functions,  such  as  power-plant  oper¬ 
ations  and  outdoor  burning  for  incineration  or  training  of  fire-fighting  crews. 

3-  Studies  of  Pollution  Potential. 

In  the  air  pollution  analysis  it  is  necessary  to  consider  the  prevailing 
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winds  and  the  modifications  that  are  caused  by  topography  or  diurnal  changes. 
Combined  with  studies  of  the  meat)  wind  speed  and  direction  should  be  studies  of 
the  diffusive  powsr  of  the  atmosphere.  In  most  oases*  the  diffuslvity  would  be 
estimated  based  on  the  considerations  of  Chapters  4  and  5.  If  a  really  critical 
determination  of  concentrations  is  required*  an  air-tracer  survey  should  be  con¬ 
ducted  at  the  site  in  order  to  provide  a  reliable  system  for  deciding  whether 
or  not  to  curtail  operations  at  a  given  time. 

The  evaluation  of  potential  hasards  both  for  site-selection  purposes  end  for 
speolflc  releases  Involves  e  knowledge  cf  atmospheric  behavior  on  all  scales 
from  the  highs  and  lows  of  the  synoptic  weather  map  to  the  individual  eddies 
that  dilute  and  redistribute  the  toxic  matter  and  other  pollutants  within  the 
atmosphere . 

Special  summaries  of  hourly  weather  observations  can  yield  a  climatology  of 
diffusion  according  to  the  Pasqulll-Oifford  system.  First*  the  frequency  of 
each  category  is  computed  for  each  wind  direction.  Next*  an  estimated  median- 
wind-speed  for  each  category  is  used  together  with  mathematical  representations 
of  Figures  21  and  22  to  arrive  at  average  concentrations. 

4.  Support  for  Current  Operations. 

Two  types  of  pollution  concentration  estimates  are  required  on  a  mo re -or-:  „ 
"real  time"  basis.  One  is  used  for  a  decision  whether  or  not  to  curtail  opera¬ 
tions.  The  other*  pertaining  only  to  toxic  substances*  is  the  area  to  evacuate 
in  case  of  an  accidental  release.  Preparing  a  radioactive  fallout  plot  is  an 
example  of  this  type  of  hazard  analysis.  Related  to  this  second  category  of 
support  is  operational  forecasting  for  intentional  releases  of  toxic  materials 
such  as  tests  of  CW  agents  or  tests  of  rocket  engines  using  toxic  propellants. 

Forecasting  for  intentional  releases  requires  more  precision  than  evaluating 
the  hazard  areas  for  accidents*  because  the  amount  of  material  released  is  known 
for  the  intentional  case,.  It  has  been  said  with  regard  to  accidental  releases 
that  "by  the  time  you  figure  out  how  much  material  was  released*  the  incident  is 
over  except  for  settling  the  lawsuito."  Therefore*  it  is  important  to  determine 
in  advance  the  probable  hazard  areas  for  a  number  of  typical  situations.  The 
situation  which  is  commonly  investigated  by  the  Atomic  Energy  Commission  is  the 
maximum  credible  accident.  A  release  of  this  magnitude  is  then  evaluated  for 
hazard  distances  under  a  variety  of  meteorological  conditions. 

5 .  Sources  of  Data  and  Technical  Assistance. 

Within  the  Air  Weather  Service,  expert  assistance  on  diffusion  and  air  pol¬ 
lution  problems  may  be  obtained  from  the  USAFETAC.  This  unit*  which  has  long 
provided  climatological  services*  haB  been  given  the  mission  of  providing  micro- 
meteorological  services  as  well.  USAFETAC  can  provide  studies  based  on  data 
from  routine  meteorological  observations  taken  almost  anywhere.  There  are  also 
limited  amounts  of  micro meteorological  data*  mostly  from  existing  missile 
ranges.  USAFETAC  personnel  can  assist  in  interpreting  these  data  as  they  apply 
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to  the  pattern  of  dispersion  of  material  released  into  the  atmosphere. 

In  certain  cases  the  site  may  be  so  remote  that  routine  meteorological  ob¬ 
servations  do  not  exist.  Also,  it  has  been  the  experience  of  many  workers  in 
micrometeorology  that  it  is  quite  difficult  to  quantitatively  predict  the  mag¬ 
nitude  of  the  local  variations  even  when  the  qualitative  variations  are  under¬ 
stood.  For  one  or  both  of  these  reasons,  it  may  prove  necessary  to  employ 
mobile  meteorological  observation  teams.  Air  Weather  Service  has  a  capability 
to  provide  routine  meteorological  observations  from  remote  sites  through  the 
resources  of  the  6th  Weather  Squadron.  If  the  requirement  for  remote  observa¬ 
tions  is  recognized  in  time  to  program  this  support  in  the  normal  budgetary 
cycle,  the  remote-site  observation  program  can  be  funded  by  Air  Weather  Service. 
If  the  requirement  is  not  recognized  in  time  for  the  normal  budgetary  processes, 
the  remote-site  observation  program  must  be  paid  for  with  a  transfer  of  funds 
to  Air  Weather  Service  by  the  operational  agency  that  needs  to  use  the  remote 
site.  Time  limitations  may  dictate  that  the  observational  program  be  handled 
by  a  contractor,  possibly  in  connection  with  an  air-tracer  survey  as  outlined 
below. 

In  cases  where  frequent  releases  of  toxic  mat>  rial  are  contemplated,  viz., 
test  ranges,  and  air-tracer  survey  should  be  undertaken  in  order  to  provide  a 
reliable  system  for  deciding  whether  or  not  to  release  at  a  given  time.  An  air- 
tracer  survey  might  even  be  Justified  to  fully  determine  the  extent  of  the  haz¬ 
ard  associated  with  the  maximum  credible  accident.  An  air-tracer  survey  con¬ 
sists  of  releasing  known  quantities  of  inert  material  into  the  atmosphere  under 
a  variety  of  meteorological  conditions  and  sampling  the  concentration  of  this 
material  at  a  variety  of  downwind  and  crosswind  distances.  The  staff  meteor¬ 
ologist  can  and  should  recommend  to  the  operator  of  the  proposed  facility  that 
an  air-tracer  survey  be  conducted  when  necessary.  The  organizations  that  con¬ 
duct  air-tracer  surveys  can  also  take  the  routine  meteorological  observations 
if  these  cannot  be  obtained  from  Air  Weather  Service  resources,  The  USAFETAC 
is  familiar  with  tracer  studies  which  have  been  conducted  in  the  past  and  can 
provide  assistance  in  defining  objectives  for  contract  meteorological  programs. 
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Appendix 

SUMMARY  OP  COMPUTATIONAL  TECHNIQUES 


In  this  Appendix  the  computational  aide  and  graphs  (full  alee)  are  assem¬ 
bled  at  the  end  for  quick  reference.  Step-by-step  procedures  are  given  for 
several  examples. 

The  systematic  nature  of  these  calculations  can  easily  give  the  impression 
that  air  pollution  concentrations  can  be  accurately  estimated.  It  should  al¬ 
ways  be  kept  in  mind  that  the  concentrations  calculated  by  these  techniques 
apply  to  open,  level  terrain.  For  an  area  characterized  by  irregular  terrain, 
the  situation  is  severely  complicated  by  the  many  local  factors  that  were  dis¬ 
cussed  in  Chapter  2.  The  estimated  concentrations  should  then  be  subjectively 
modified  to  account  for  local  changes  in  the  flow  pattern  and  differences  in 
the  Intensity  of  turbulent  mixing. 

1.  General  Form  of  Calculations. 

The  concentration  of  pollutants  consists  of  the  background  level  plus  the 
contribution  due  to  local  sources.  All  of  the  methods  discussed  in  Chapters  4 
and  5  concern  diffusion  from  local  sources.  For  simplicity,  the  background 
concentration  will  be  omitted  in  the  examples  Which  follow. 

The  contribution  due  to  local  sources  is  computed  by  evaluating  certain 
factors  in  order.  Each  additional  consideration  can  be  accounted  for  by  multi¬ 
plying  the  result  of  the  previous  step  by  another  factor.  The  complete  equa¬ 
tion  for  concentration  from  a  single  source  can  be  represented  as  a  series  of 
products,  thus 

q  ■  [Source  Strength]  x  [Along-Wind  Dilution]  x  [Cross-Section  Dilution] 
x  [Vertical  Term]  x  [Lateral  Term]  x  [Depletion  Term]. 

Any  time  a  factor  can  be  neglected,  the  value  of  its  term  is  unity.  Since  mul¬ 
tiplying  by  one  leaves  the  answer  unchanged,  the  unneeded  terms  are  simply 
omitted.  The  first  three  terms  must  always  be  used,  and  the  vertical  term  is 
usually  needed.  Numerical  values  of  the  various  factors  can  be  obtained  from 
Figures  21  through  34,  36,  37,  47,  and  48  with  use  of  Tables  1  and  2.  The  use 
of  these  values  in  the  equation  will  now  be  illustrated. 

a.  Source  Strength.  If  the  source  is  continuous .  it  should  be  stated  in 
grams  per  second.  If  the  source  is  instantaneous,  the  source  strength  is  simply 
the  mass  of  contaminant  release.  All  of  the  graphical  aids  in  this  report  re¬ 
quire  metric  units,  i.e.,  meters,  seconds,  grams,  etc.  Thus,  the  mass  should 
be  converted  to  grams  (one  pound  equals  453.6  grams). 
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b.  Alona-Wind  Dilution. 

(1)  Continuous  Source.  Divide  the  source  strength  by  the  wind  epeed 
(meters  per  seoond)  (1  kt  ■  0.515  m/sec).  For  smokestacks,  the  wind  at  release 
(stack)  height  can  be  estimated  from  Figure  49.  For  example,  the  wind  at  66 
feet  would  be  1.33  times  the  strength  of  the  anemometer-level  wind  under  neu¬ 
tral  conditions. 


Figure  49.  Graph  for  Determining  Wind  Speed  at  Release 
(stack)  Height. 

(2)  Instantaneous  Source.  Select  a  value  from  the  right-hand  scale  of 
Figure  32  according  to  the  applicable  stability  category  and  the  desired  dis¬ 
tance  downwind  of  the  release  point.  For  example,  at  1000  meters  downwind  in 
neutral  conditions,  the  divisor  would  be  about  280  on  the  logarithmic  scale. 
Select  the  appropriate  line  as  follows!  If  the  stability  category  from  Table  1 
and  Table  2  is  A  or  B,  use  the  line  marked  "unstable."  The  "neutral"  line 
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corresponds  to  Category  D  and  the  "stable"  line  is  equivalent  to  Category  E  or 
F.  For  Category  C,  use  a  value  between  the  "unstable"  and  "neutral"  values. 

c.  Cross-Section  Dilution.  The  vertieal  and  oross-wind  dimensions  of  the 
plume  or  puff  increase  downwind.  The  same  amount  of  material  must  occupy  a 
larger  and  larger  space,  thus  lowering  the  concentration  at  the  center  of  the 
plume  or  puff.  After  this  step  the  units  will  be  grams  per  oubio  meter.  Pro¬ 
ceed  as  follows t 

(1)  Continuous  Source.  Use  the  value  on  the  ordinate  of  Figure  23  as 
a  multiplier. 

(2)  Instantaneous  Source.  Use  the  value  on  the  ordinate  of  Figure  34 
as  a  multiplier.  For  example,  at  a  downwind  distance  of  300  meters  in  unstable 
conditions,  the  factor  would  be  about  5  x  10"*.  Selection  of  the  three  cate¬ 
gories  waB  explained  in  paragraph  1(b)(2)  of  this  Appendix. 

d.  The  Vertical  Term.  The  vertical  term  is  needed  in  most  cases,  even  for 
concentrations  measured  at  the  height  of  the  cloud  center.  In  its  most  general 
form,  the  vertical  term  is 

-(z  -  h)’/2o  *  -(z  +  h)*/2o  • 

e  +  e 

where  h  is  the  height  of  the  cloud  center  above  ground  arid  z  is  the  height  above 
ground  where  the  concentration  is  measured.  This  general  form  is  given  here 
only  for  the  sake  of  completeness  and  is  not  used  in  the  examples  shown  in  this 
Appendix.  An  explanation  of  the  equation  and  the  reason  for  the  second  term  is 
given  by  Sutton  ([22]  pp.  139-l4o)  and  in  Meteorology  and  Atomic  Energy  1968 

([20]  p.  99). 

When  measuring  at  ground  level,  z  becomes  zero  and  the  general  equation  re¬ 
duces  to 

-ha/2o  8 
2e  z 

with  h,  the  height  of  the  cloud  center  above  the  ground.  It  is  this  reduced 
form  of  the  equation  that  is  used  in  the  examples. 

To  aid  in  the  calculations.  Figure  24  was  developed.  It  gives  values  of 
-a8 /2ba 

e ■  '  on  the  ordinate  for  any  value  of  a  on  the  abscissa.  Values  of  the 

parameter  b  are  shown  as  a  family  of  sloping  lines.  To  use  th&  figures  to 
evaluate  the  vertical  term,  use  the  value  of  h  as  the  abscissa  and  proceed  up¬ 
ward  to  the  (interpolated)  sloping  line  corresponding  to  &'z  and  proceed  across 
to  the  appropriate  value  of 

-ha/2o  8 
e 

on  the  ordinate.  Tnen  double  the  result. 
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The  vertioel  term  must  bo  modified  if  buoyant  aouroes  are  involved*  or  if 
there  la  a  aignlfioant  elevated  inveralon.  If  the  aouroe  la  buoyant*  the  ef¬ 
fective  eouroe  height  muet  be  found.  This  modified  value  of  h  ia  then  uaed  to 
evaluate  the  vertical  term  in  the  uaual  way.  Rulea-of -thumb  for  eatimating  the 
effect  of  inveralona  are  given  in  Chapter  4.  The  neceaaary  atepa  for  varloua 
aituationa  can  be  aummarlaed  aa  follows t 

(1)  Non-Buoyant  Sources  without  Elevated  Inversions.  For  a  given  sta¬ 
bility  category  and  diatanoe  downwind,  find  o*  from  Figure  22.  Evaluate  the 


term 


2e 


-ha/2 o  ■ 

Z 


with  the  aid  of  Figure  24*  aa  described  above. 

(2)  Buoyant  Sources  without  Elevated  Inveralona.  Calculate  Ah  accord¬ 
ing  to  the  Holland  formula 


Ah 


3  V  d  T 

— tjS—  [l.5  +  2.68  x  10"a  n  (l  -  £t)  d] 

s 


Next,  compute  the  buoyancy  flux  F  «  g  Vs  Ra  (1  -  Ta/Tfl).  (The  definitions  of 
the  symbols  were  given  in  Chapter  5.)  Again*  all  units  should  be  in  meters  or 
meters  per  second.  Enter  Figure  47  with  the  buoyancy  flux  as  the  abscissa  and 
proceed  up  to  the  sloping  line  corresponding  to  the  wind  speed.  From  this 
point,  read  across  to  the  ordinate  to  find  the  value  2F^*  u"1 .  Now,  enter  Fig¬ 
ure  48  with  downwind  distance  as  the  abscissa  and  proceed  up  to  the  sloping 
line  corresponding  to  2F^*  U  1 .  From  that  point*  proceed  across  to  the  ordinate 
to  find  a  second  value  of  Ah.  Use  whichever  value  of  Ah  is  smaller.  Add  Ah  to 
h  and  then  proceed  as  in  (1)  above. 

(3)  Significant  Inversion  Aloft.  First,  use  Figure  22  to  find  oz. 

There  are  three  types  of  situations j 

(a)  For  a>z  less  than  half  the  inversion  height,  neglect  the  in¬ 
version  and  proceed  as  in  (1)  above. 

(b)  For  oz  more  than  twice  the  inversion  height,  the  vertical 
term  is  2.56  crz/H,  where  H  i3  the  inversion  height.  By  using  this  vertical 
term,  we  convert  Equation  (10)  to  Equation  (17). 

(c)  For  oz  between  half  and  twice  the  inversion  height,  proceed 
as  follows.  First  evaluate  the  vertical  term  according  to  (1)  or  (2)  above  for 
az  equal  to  half  the  Inversion  height.  Plot  this  value  on  a  log-log  plot  (such 
as  Figure  22  or  23)  at  the  downwind  distance  where  cr^  is  half  the  inversion 
height.  At  the  distance  where  is  twice  the  inversion  height,  plot  the  value 
5.12.  The  vertical  term  fox  any  downwind  distance  in  this  range  can  be  found 
by  Joining  the  two  plotted  points  with  a  straight  line.  These  steps  are  essen¬ 
tially  equivalent  to  following  the  "rules  of  thumb"  given  in  Chapter  4,  para¬ 


graph  7. 


There  would  be  a  slight  difference  because  the  plot  of  o  on  log-log 

z 
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paper  is  not  a  atraight  line. 

e.  The  Lateral  Term.  Thia  term  la  not  needed  for  calculating  peak  oon- 
centratlona  unless  the  peak  is  the  result  of  the  combined  effect  of  several 
scattered  sources.  The  value  of  the  lateral  term  is  always  unity  for  peak  con¬ 
centrations  from  single  sources.  The  form  of  the  term  is 

-y"/2c  * 
e  y 

where  y  is  the  cross-wind  distance  in  meters. 

(1)  Continuous  Sources.  Figure  21  gives  values  of  oy  as  a  function  of 
downwind  distances  and  stability  category  for  the  Pasquill-Qifford  System.  Use 
Figure  21  to  get  0 and  then  use  Figure  24  to  evaluate 

-ya/2o  * 
e  * 

as  explained  for  the  vertical  term.  Isopleths  of  q  U/Q  for  ground-level  concen¬ 
trations  from  ground-level  sources  are  provided  as  Figure  25  through  30.  The 
concentration,  off  axis,  for  any  other  values  of  z  and  h  can  be  found  by  using 
half  the  vertical  term  as  a  multiplier  for  the  values  in  Figures  25  through  30. 

(2)  Instantaneous  Sources.  Figure  33  gives  values  of  cty  for  three 
stability  categories:  unstable,  neutral,  and  stable.  Select  the  stobility 
category  from  Table  2  and  the  equivalence  suggested  for  "along-wlnd  dilution." 
Use  these  values  of  oy  with  Figure  24  to  evaluate  the  lateral  term. 

f •  The  Depletion  Term.  Various  depletion  processes  are  discussed  at  the 
beginning  of  Chapter  5. 

(1)  Concentrations  Remaining  in  Pollutant  Cloud.  Calculate  depletion 
rate  per  meter  (K/U  where  K  is  depletion  per  unit  time).  Be  sure  the  same  time 
units  are  used  for  U  and  K.  Wind  speeds  in  knots  can  bs  used  with  Figure  3 6  to 
find  K/U  if  K  is  expressed  in  depletion  per  minute.  Figure  37  gives  the  deple¬ 
tion  term  as  a  func  '.on  of  distance  for  a  family  of  K/U  valueB. 

(2)  Pollutants  Deposited  on  the  Ground.  There  are  two  general  cases: 
dry  deposition  and  washout  due  to  rain.  In  the  dry  case,  the  deposit  rate  is 
proportional  to  the  ground-level  concentration.  For  washout,  the  deposit  is 
proportional  to  the  vertically-integrated  depletion.  Refer  to  Chapter  5  for 
details . 

2.  Illustrative  Examples. 

In  this  section  we  present  some  complicated  problems  which  require  evalua¬ 
tion  of  all  terms  except  the  depletion  term  and  show  the  step-by-step  solution 
of  these  problems. 

Example  1.  It  is  a  summer  afternoon  at  Edwards  AFB  and  the  surface  wind  is 
from  350°  at  13  kts.  A  small  experimental  rocket  motor  containing  20  grams  of 
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finely  divided  beryllium  explodes  on  the  top  of  a  15-ft  test  stand  releasing  a 
cloud  of  beryllium  particles  which,  for  this  problem,  Is  assumed  to  have  no 
buoyancy.  A  group  of  engineers  was  standing  150  yards  south  of  the  test  stand 
at  the  time  of  the  explosion.  What  was  the  maximum  concentration  of  beryllium 
that  the  engineers  were  exposed  to? 


Figure  50.  Diagram  of  Example  1. 


Solution: 

Step  1.  Source  strength:  This  is  given  as  20  grams j  conversion  to  other  units 
is  not  required  in  this  case. 

Step  2.  Along -wind  dilution  factor:  Using  Figure  32  find  the  horizontal  down¬ 
wind  distance,  x  (in  meters). 

x  ■  150  yds  x  3'.g^  ft^meter  x  coslne  (10° )  "  135  meters 

Using  this  distance,  enter  Figure  32  to  get  the  downwind  dilution  fac¬ 
tor.  Table  2  gives  the  stability  category  as  C.  Interpolate  between 
the  "neutral"  and  "unstable"  line  to  obtain  25  meters  as  the  factor. 
Dividing  20  by  25  gives  0.80  grams/meter. 

Step  3.  Cross-section  dilution:  Go  to  Figure  34.  For  a  downwind  distance  of 
135  meters,  interpolate  between  the  "neutral"  and  "unstable"  curves  to 
obtain  the  factor  4.6  x  10"S  per  square  meter.  Multiply  this  by  the 
result  of  Step  2  to  obtain  4.6  x  I0“a  x  0.30  «  3.68  x  I0”a  grams  per 
cubic  meter  on  the  centerline. 

Step  4.  Vertical  term:  The  Pasquill  stability  category  C  has  already  been 
selected  from  Table  2.  For  this  category.  Figure  22  gives  oz  *  10m 
for  the  downwind  distance  135  meters.  There  is  no  limiting  inversion 
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and  we  are  neglecting  buoyancy,  thus,  the  vertical  term  need  not  be  . 
modified.  Convert  h  to  meters,  i.e.,  h  ■  15/3.28  »  4.57  m.  Using 
Figure  24  with  a  *  4.57  and  b  »  10,  the  vertical  term  is  2  x  0.90  « 
1.80.  Thus,  the  concentration  directly  under  the  cloud  would  be 
1.80  x  3.68  x  10" 3  -  6.62  x  io“a  g/ma . 

Step  5.  Lateral  term:  The  lateral  distance  from  the  cloud  trajectory  is 

150  yd  x  3  yd/3.28  m  x  Sin  (10°)  «  24  meters.  Enter  Figure  33  with 
the  downwind  distance  (135  m)  and  interpolate  again  between  the  "neu¬ 
tral"  and  "unstable"  lines  to  obtain  oy,  which  is  7.5  meters.  Use 
Figure  24  with  a  *»  24  and  b  «*  7-5  to  get  the  lateral  term,  which  is 
0.006.  Multiplying  this  by  the  result  of  Step  4  gives  the  answer, 

q"4  x  io'B  grams/ma 

Example  2.  The  base  civil  engineer  proposes  that  a  new  base  heating  plant  will 
have  a  65-ft  stack,  four  feet  in  diameter  at  the  top.  The  stack  temperature 
will  be  600° F  and  stack  velocity  10  ft/sec.  Sulfur  in  the  fuel  oil  is  expected 
to  produce  270  grams/sec  of  S08  as  a  combustion  product.  If  the  dining  hall  is 
2000  feet  northeast  of  the  proposed  heating  plant,  what  would  be  the  concentra¬ 
tion  of  S08  on  the  ground  at  the  dining  hall  at  sunrise  on  a  clear  morning  if 
the  winds  were  southwest  at  4  kts  at  anemometer  level?  Assume  the  air  tempera¬ 
ture  at  stack  height  is  50°F  and  the  station  pressure  is  1000  mb. 


Figure  51.  Diagram  of  Example  2. 
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Solution; 

Step  1.  Convert  all  units  to  metric.  The  stack  height  is  65/3.28  »  20  meters, 
the  stack  diameter  ■  1.22  meters,  the  stack  velocity  3  m/aect  the 
stack  temperature  3l6°C  ■  589° K,  while  the  air  temperature  is  283°K. 
The  surface  (anemometer)  wind  is  2.06  m/sec  and  the  downwind  distance 
to  the  dining  hall  is  610  meters.  The  continuous  source  strength  (2?0 
gm/sec)  is  already  in  metric  units,  and  the  stability  condition  is 
stable  (Category  E). 

Step  2.  Along-wind  dilution;  First,  find  the  wind  at  stack  height.,  From 

Figure  49,  the  multiplier  would  be  I.95  under  stable  conditions,  thus, 
the  wind  speed  is  4  m/sec.  Dividing  the  release  rata  (270  grams/sec) 
by  4  m/sec  gives  67.5  grams/meter. 

Step  3.  Cross-section  dilution;  The  stability  category  determined  from  Tables 
1  and  2  is  E.  The  dilution  factor.  Figure  23,  for  Category  E  and  a 
downwind  distance  of  610  meters  is  3-43  x  10“"*.  Multiply  this  by  67.5 
(the  result  of  Step  2)  to  get  2.30  x  10"3  g/m8 . 

Step  4.  Vertical  term:  This  problem  involves  a  buoyant  source,  so  the  first 
thing  is  to  evaluate  the  additional  rise,  Ah,  above  the  stack  height. 
The  Holland  formula  for  Ah  with  a  correction  of  3  gives 

Ah  =  [1.5  +  2.68  x  10  “*  x  1000  ^1.0  -  ||$)(1.&2)] 

-  2.74  (1.5  +  1.71)  -  8.8  m 

Compare  this  with  the  "transitional"  plume  rise  given  by  Figures  47 
and  48.  The  buoyancy  flux,  F,  is 

9.0  x  3.0  *  (.61)*  (l.o  -fjg)  -  5.T 

Enter  Figure  47  with  5*7  on  the  abscissa  and  read  up  to  the  4  n\/sec 
(stack-height  wind)  line.  The  value  from  thiB  graph  (read  from  the 
ordinate)  is  0.9,  Enter  Figure  48  et  a  downwind  distance  of  610  meters 
and  proceed  up  to  the  (interpolated)  sloping  line  corresponding  to  0.9 
from  Figure  47.  Reading  across  to  the  ordinate,  find  Ah  «  64  meters. 
Use  the  lesser  value  for  Ah.  Add  this  to  20  m  (the  original  stack 
height)  to  give  an  effective  height  of  29  meters.  Next,  find  cz  from 
Figure  22.  In  this  case,  c£  is  14  meters.  At  ground  level  2  •»  0. 

ThUB,  |z+h|=  |z  -  h|  -  29.  The  vertical  term  (using  Figure  24  with 
a  =  29  and  b  »  14)  becomes  0.1  +  0.1  =  0.2,  Thus,  the  ground-level 
concentration  is  4.6  x  I0“a  grams/m8  or  4.6  milligrams  per  cubic  meter. 

Example  3.  Suppose  the  heating  plant  of  Example  2  continued  to  emit  270  grams/ 
sec  of  S0»  two  hours  after  sunrise,  at  which  time  the  ground-based  inversion 
had  lifted  to  200  feet  and  the  surface  winds  had  increased  to  9  knots.  What 
would  the  S0B  concentration  at  the  dining  hall  be  then? 
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Figure  52.  Diagram  of  Example  3. 

Solution: 

Step  1.  All  values  are  the  same  as  it:  the  last  problem  except  the  surface  wind 
is  4.6  m/sec  and  the  stability  category  is  now  C. 

Step  2.  The  rate  of  increase  of  wind  with  height  should  now  be  interpolated 

between  the  "neutral"  and  "unstable"  lines  of  Figure  49,  giving  us  now 
an  increase  factor  of  1.3>  giving  6.0  nysec.  Dividing  270  by  6.0  gives 
45.0  grams/meter. 

Step  3.  The  cross-section  dilution  factor  for  610  meters  and  C  stability  is 
6.0  x  io“\  Thus,  the  concentration  on  the  centerline  of  the  smoke 
plume  is  45.0  x  6.0  x  10-8  =  2,7  x  10-3  grams/m*. 

Step  4.  Find  the  value  of ' og  corresponding  to  6l0  meters  downwind  for  C  stabil¬ 
ity.  From  Figure  22,  is  37  meters.  This  is  mors  than  half  the  in¬ 
version  height  (61  meters),  but  less  than  twice  the  inversion  height. 
Thus,  the  method  given  in  paragraph  ld(3)  of  the  Appendix  must  be  used. 
The  distance  at  v/hich  a would  be  half  the  inversion  height  is  48o 
meters  downwind  for  category  C  stability.  Similarly.-  Figure  22  gives 
2200  meters  as  the  distance  where  c  would  be  twice  the  inversion 

Z 

height.  First,  evaluate  the  vertical  term  in  the  normal  way  for  e  o„ 

z 

of  30  meters.  The  effective  source  height  is  now  somewhat  lower  be¬ 
cause  of  the  increased  wind  speed  at  stack  height.  Ah  is  now 

(1.50  +  1.71)  «  5.87  m 
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This;  added  to  the  20-mater  stack  height,  makes  the  affect 4  ve  height 
now  26  metefs.  Enter  Figure  24  with  a  ■  2t3  and  b  ■  30/  getting  0.72 
for  the  exponential  term,  This  term,  in  turn,  if  doubled  to  get  the 
vortical  term  (■  1;44)  at  430  meters  iownw/vrt.  Next,  find  any  conven¬ 
ient  log-log  paper.  (Figure  21  or  Figure  22  would  do.)  .'At  480  meters 
downwind,  plot  k  point  with  an  ordinate,  value  of  1.44.  At  2.9  km  down¬ 
wind,  plot  a  point  with  5 .12  as  the  ordinate .  Join  thede  two  points 
with  a  straight  line  and!  read  l.?5  as  the  vertical  term  for  a  downwind 
distance  of  6iO  meters.  Multiply  1.75  by  the  result  of  Step  3s 

S'  ,  .  ’  '<  ■  1  .  t  ■  . 

2.7  x  10“*  x  1.75 >  4.7  x  10“a  grams/m® 

Notice  that  this  concant ration  it  higher  than  the  concentration  calcu¬ 
lated  in  Example  2,  even , though  the  surface  wind  is  stronger  and  the 
lapse  rate  is  unstable  next  ,tb  the  ground. 


3. 


Computational  Tables  andQraphs.  ,  , 

Tables  1  and  2  along  with  thi  figures  appearing  on  the  following  pages  are 
repeated  from  preceding  chapters  of  this  .guide  to  allow  easy  access  to  the  in¬ 
formation  required  for  the  calculations  within  this  Appendix.  The  larger  size ■. 
of  these  duplicated  graphs  will  allow  more  accurate  data  to  be  read  and  their 
use  in  solving  problems  in  the  field  will  be  enhanced. 

’  '  /•  l  '  ,  ‘  • 

* '  •  , 

TABLE  1 

Net  Radiation  Classes 


_ „ _ 

Elevation  of  Sun 

Cloud  Cover 

Night 

*  15° 

>  15°  but 
*  35° 

>  35°  but 
<  65° 

>  65° 

<  3/8 

-2 

-1 

+1 

+2 

+3 

3/8  or  4/8  at  any 
hgt  or  broken  abv 
16000  ft  (incl.  -©) 

-1 

0 

+1 

+2 

+3 

©  abv  16000  ft  or 
®  between  7000  ft 
and  16000  ft 

-1 

0 

0 

+1 

+1 

®  below  7000  ft 

0 

0 

0 

0 

+1 

©  below  7000  ft 

0 

0 

0 

0 

0 

gk 
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TABLE  2 

Stability  Category  Baaed  On  Wind  Speed  and  Net  Radiation 
Surface  Wind  Speed  Net  Radiation  Olaee  (see  Table  1) 
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